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ABSTRACT 


A surface-condenser analysis coda was coupled with a 
constrained function minimization soda to produce an auto- 
mated marine condenser design and optimization package. The 
program, CONDIP, was based on tha principles developed in 
ORCON1, a sophisticated computer soda produced by «he Oak 
Ridge National Laboratory. CONMIN, the optimizacion progzraa, 
was developed at the Ames Research Center. 

CCNDIP is an extremely versatile design tool, incerpo- 
rating a detailed analysis of the complex steam-side thermo- 
dynamic processes occurring at each row in the condenser. 
The additicnal capability of tube enhancement is also 
included. However, in coupling CONDIP with CONMIN numerous 
problems had to be overcom? in ordar to make CONDIP capable 
of ccapleting an analysis even whan theracdynamic conditicns 
in the condenser became infeasible. This ha@ ‘+o be accon- 
plished whila ensuring continuity in ail constraint and 
objective furction evaluations. A s2rias of *est cas3s were 
conducted to evaluate ani compare tha importance of various 
objective functions and design critaria. 
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I. [NPRODUCTION 


A. BACKROUND 


Por aany years the steam plant was unchallenged in its 
tole as the primary type of marina propulsion systems. But 
recently gas turbines have become a desirable alternative 
despite the fact that they are lass sfficient ‘han compa- 
rable steam plants. The primary aivantage of gas turbines 
is their light weight and compact size. Thus, in crder for 
the marine steam plant to survive, it is imperative that 

aghter, more compact and efficient steam plaats be 
developed. 

While there are numerous advanced concepts in all areas 
of steam prepulsion which can be explored, one simple way to 
streamline the steam plant is the elimination of overdesign. 
Most overdesign is due to unnecessary safety factors used to 
offset lack of detailed knowledge about the thermai 
processes in the plant. I@entification of the ainimun safe 
design could significantly reduce unnecessary overdesign and 
result in the development of a smaller, more compact power 
plant. 


B. HSETHODOLOGY 


In the United States the most prevalent criterion for 
the design and specification of surface condensers is based 
on the "square-root of ¥" relationship, as developed by the 
Heat Exhange Institute (HEI) (Ref. 1]. The HEI method was 
adopted by the Department of the Navy Bureau of Ships (now 
Naval Sea Systern Coaaand) for specification of U.S. Naval 
condensers. 
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The HEI method is very simple 1n its appreach, calcu- 
lating tke overall heat transfer coefficient as a funccior 
of cooling water vslocity through the tubes, inlet coclant 
temprature, tube wall thickness ani material, and fouling. 
The limitations of this method ar3 apparent. Designs based 
on HEI are insensitive to shell-sida conditions. Saturation 

team pressure, temperature and enthalpy are assumed +o be 

constant as steam passes through the bundle, whereas in 
reality there is a continual pressure drop as steam flow 
passes over the rows of tubes, with a corresponding decrease 
in saturaticn temperature. There is also no previsior for 
any effects of condensate film, external tube enhancement, 
etc. on the shell-side of the bundle. In addition, the HEI 
method does not account for the prasence and effect of non- 
condensible gases that inevitably contaminate a condenser. 

With the capabilities of high speed computers now avail- 
abl3, more comprehensive methods have been developed to 
account for the daficiencies of the HEI method. In pazrtic- 
ular, a radial flow computer code was developed to calculate 
the local heat transfer and thermodynamic properties ona 
row by row basis. Known as ORCON1, chis code was developed 
by Oak Ridge National Laboratory (ORNL) under contract +o 
the Office cf Saline Water during th2 period from 1968-1970 
(Ref. 2). The program was based, in part, on the work 
performed by Zissenberg (Ref. 3}. Eissenberg's experimenzal 
results led to correction factors on the basic Nusselt equa- 
tion to acccunt for condensate inundation effects on tubes 
within a condenser bundle. Basically, ORCON1 divides the 
condenser into sectors and performs a row by row analysis 
within each sector, determing local heat ‘¢ransfer coeffi- 
cients, heat flux, steam charactaristics, the affect of 
condensate inundation and numerous other parameters at each 
row. ORCON1 is also capable of incorporating the effects of 
both tube-side and steam-side enhancement factors. Since 
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ORCON1 represented a auch sore comprehensive and detailed 
analysis of -the condenser then tha less exact HEI nethod, 
its results could be expacted to be more precise. 

Some work has been done at the Naval Postgraduate School 
to improve +he capabilities of ORCON1. In his ievelopement 
of OPCODE2, Johnson (Ref. 5] added subroutines <0 ORCON1 
which calculated tube-side ressure drops, corresponding 
pumping power and condenser volume. Nunn and Marto [Ref. 14] 
have incorporated the effects of vapor shear in ar amended 
version of ORCON1 called MORCON. MORCON includes the corre- 
lations developed by Fujii [Ref. 4] to determine tke effec 
of vapor velocity on the thermal resistance of the conden- 
sate film on the condenser tubes. In general, vapor shear 
effects tend to enhance the condenser. heat transfer on the 
steam-side of the tube while condensate inundation tends to 
inhibit it. 

The ability to represent numerically the actual thermo- 
dynamic processes occurring within the condenser has 
improved dramatically. However, the capability to couple 
these increasingly comprehensive and complex conderser 
design programs with an optimizing procedure has not made 
comparable progress. Optimization is a powerful tool which 
can help in reducing overdesign and achieving the goal of a 
safe compact condenser design. 

There are currently numerous computer optimization 
programs available which can be coupled with general design 
progrags of all types to numerically improve and ultimately 
determine the best design. The key is to properly write the 
design program so that it is compatible with the optiaizer. 
Johnson (Ref. 5] developed a computer program called 
OPCODE1, based on the HEI sethod of condenser design, and 
was able to couple it with one such numerical optimizer. The 
results of OPCODE1 demonstrated how condenser designs can 
indeed be safely improved upon. It also revealad the 
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wersatility of condenser design optimization as a powerful 
design tool. However, Johnson was unsuccessful in coupling 
OPCODE2 (his derivative of ORCON1) with an optimizer. This 
failure does not alter the fact that that in order to fully 
appreciate more sophisticated condenser design analyses, 
such as that used in ORCON1, it is imperative that coaputer 
programs be developed which will be compatible with current 
humerical optimizers. 


C. OBJECTIVE 


There were two primary objectives of this thesis. The 
first objective was to develop a computer code which incor- 
porates the basic ccndenser analysis of ORCON1 and the 
subsequent iaprovements made in MIRCON and OPCODE2, but 
which will be capable of being coupled with a numerical 
optimizer to yield a complete, detailed design package. 
This design package can then be usad as a tool in obtaining 
a much sore reasonable conceptual dasign and for use in 
comparison studies. It would provide the naval architect 
the ability to optimize weight, volume, cost or any other 
potential design objective of the aarine plant. . 

The second objective was to make this design package 
capable of determining the single best design rather than 
simply an improvement over the initial design. The key was 
to construct the program in such a way so the optiaizer does 
not stop at some relative optimus, but continues the anal- 
ysis until no further iaprovement scan be realized. It is 
most desirable to be able to reach this single true optiniun 
design regardless of initial design variable values. 
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TI. NUSEBICAL OPTINIZATION 


A. BACKGROUND 


Nearly all design problems require either the miniaiza- 
tion or maximization of a parameter. This parameter will be 
called che problem's cbhjective function or design objective 
(Ref. 6}. Por a given design to ba feasible or acceptable, 
it aust satisfy a set of design constraints which are either 
maximum or minimum lisgiting values for a pra-deterzined set 
of parameters or functions of parazeters. For exaaple, in 
any condenser design the outer diameter of a condenser tube 
can never be less than zero and there is normally some prac- 
tical upper limit which also cannot be exceeded. These 
limits.are design constraints on tha tube outer diameter. In 
the design problem there is also a sat of design variables 
which are parameters whose values can be changed within 
specified limits in order to minimize or maximiza the design 
objective. For exagaple, in minisizing the condenser volume 
an engineer may want to vary tube inner diametar, tube wall 
thickness and tube length. These three parameters would thus 
be exaaples of typical design variables. 

For such complex design probleas as the treatment of the 
condenser design in ORCON1, it is necessary to choose an 
optimization scheme which can handle the problerm and provide 
a rational, rapid approach to design automation and opti- 
mization. An optimization program based on direct methods 
for scluticn of constrained problaas (Ref. 11] was chosen 
for this research work. 
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B. CONSTRAINED FUNCTION HIWINIZATION (CONMIN) 


Vanderplaats (Ref. 7] developed an optimization prograa, 
CONMIN, capable of optimizing a very wide class of engi- 
neering problems. CONAIN is a fortran program, in subprogran 
ee fora, that optisizes a sulti-wariable function subject toa 
set of inequality constraints. 

It is practical at this point to introduce three basic 
definitions and their respective conditions (Ref. 8]. 


DESIGN "Het Those rametars which the optiagzacior 
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; jescciee that the optimization process requires tha 
, ainigization of a particular objective function the general 
optiazation problem can be stated as: 


torgiatsize fue sojecti¥e Fonction, Pe 
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G(x) $ 0.0 j= 1,NCON (eqn 2.1) 


VL8j 3X S VOB; i= 1,0DV (eqn 2.2) 


In the general probles, 65 (0) are the constraint 
functions; there are MNCON constraints and wWDV design 
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variables; VLB; and VUB, are the lower bounds and upper bcunds 
of the i-th design variable. If the equality conditicrn is 
met, (Gj (Z)0.), the constraint is active. If the inequalit 

is net, (G 5() <0-), the constraint is inactive. FPinally, if 
the inequality of equation 2.1 is violated, (G 5(K) 20.) that 
constraint is said to be violated. Because of numerical 
inaccuracies representing exact zero on the cocaputer, the 
equality condition is represented by a bard around the value 
G 5 (x) #(0.2C7) where CT is the constraint thickness. 

Any design which satisfies the inequalities sf equations 
2.1 and 2.2, thus hawing no violated constraints, is said to 
te feasible. If the design violates any of these constraints 
it is said to be an infeasible design. The design which best 
minimizes the objective function whila still remaining 
feasible is said to be optinal. 

; CONSIN requires an initial set of values for “he design 
wariablas xX to obtain an initial design which is either 
feasible or infeasible. If the initial design is feasible, 
CONMIN mowes in a direction which will sminimize the objec- 
tive function. If the initial design is infeasible, CONMIN 
moves toward a feasible solution with aininaal increase in 
the object function. 

The optimization process proceeds in an iterative 
fashion. Johnson [Ref. 5} presents in greater detail the 
procedures utilized in CONMIN to search for the sininaua 
objective value. In general, the methods used by CONSIN to 
Getermine search direction include the sethod of steepest 
descent, the method of conjugate direction, and the method 
of feasible directions. For further background concerning 
CONMNIN and the numerical techniques utilized in optimiza- 
tion, consult Vanderplaats (Ref. 7], Fletcher and Reeves 
(Ref. 9}, Zoutendijk (Ref. 10], and Vanderplaats and Hoses 
(ref. 12}. However, it is necessary to stress a few 
pertinent poiats which will aid in understanding how the 
program was developed in this thesis. 
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The optimization process begins by calculating the 
gradient of the objective function using a finite difference 
technique. A perturbation is appliad +to each of the design 
watiables in a single forward step and the gradient vector 
is deternined. 


p< 28») 
OF (x) = 
agsn 
The search direction is then calculated and is a function of 
this gradient and any active or violated constraints 
resulting from the applied perturbation. Subsequent search 
directions are a function of previous search directions, as 
well as current gradient information and any appropriate 
constraint factors. Obviously, the size of the perturbation 
and the size of the bandwidth about an active constraint 
will have a great deal of effect on the search direction and 
ultimat2 optimization process. This detail will be recalled 
later-on during the code development. 
‘% There are some limitations to CONMIN. The nuagaber of 
" “4 design variables (NDV) directly affects the computational 
‘| time required to reach an optimum. Since +he calculation of 
i the gradient inforgation required for each design variable 


at the beginning of each design itaration is found by using 
a single forward finite difference step, requiring a 
| | complete pass through the analysis portion of the prograa, 
i there is an increase in CPO time as NWDV inceases. Also, as 
A WOW increases, there is the corresponding rise in machine 
; related nuserical innacuracy. Vanderplaats [Ref. 6] recoa- 
" sends no more than tventy as a practical limit for the 
number of design variables. 
: It is quite possible that while design isprovement say 
i be obtained, the single best design optiauam or true optinus 
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may not te reached. This is not an uncogmon occurrence and 
there are several possible explanations. For exaaple, the 
@esign proklem may net be forsulated properly or the anal- 
ysis may be extremely complex and non-linear. However, 4 
more ccmmon reason is that there are “relative optiauas" 
between the initial design and the single true optisaua. This 
concept cf relative vs true optigauma design can be better 
explained through an analogy. The search for the best 
optimuam design can be likened to a blind man clirnbing to the 
top of a macuntain. The blind man knows he is proceeding up 
the mountain by sensing the direction of ascent. However, 
the paths he takes say be lisitai by barriers or fences 
which will restrict the directions he can go. These fences 
represent constraints in the optiaization problem. During 
the journey he may also encounter saall crests and valleys. 
If the available paths lead the blind aan up to one of these 
crests pricr to reaching th? mountain top, he will be 
confronted with a situation where he will sersa no further 
rate of ascent and he will stop his journey. So although he 
has made progress from his initial starting point, the san 
did not achieve his ultimate goal of climbing tha aountain. 
During optimization, the search for a true optimus say 
proceed along a path cn which the objective function assumes 
such relative optimus values. If the optiaizer can not be 
gpade to “look beyond" these crelative peaks, then the opti- 
mization will cease - at a design which may be an iaprove- 
gent over the initial one but short of the true optinua. 
This protlem say be overcome by starting the design with 
several different initial design vectors, X, until the sare 
optiaal design is repeated. <Anothar alternative say be to 
increase the size of the finite difference so that the opti- 
mizer uses larger perturbations of X thus looking beyond any 
small increases in the objective function which could stand 
in the way of further design progress. This second 
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alternative will be specifically addressed during the 
discussion of the code developaent. 


C. CONTROL PROGRAN POR ENGINEERING SYNTHESIS (COPES) 


The optisizer, CONMIN, was written in subroutine fora. 
Vanderplaats [Ref. 13] has developed a main program to 
simplify the use of CONMIN and aid in the design optiniza- 
tion process. 

The user sust supply an analysis subroutine called 
ANALIZ, which consists of thrae segaents: input, analysis 
and output. COPES acts as an interface between ANALIZ and 
the optimizer CONSHIN. Based on a flag from COPES 
(ICALC#1,2,3) AWALIZ performs the proper function. Figure 
2.1 offers a simplified illustration of the interrelation- 
ship bLetween COPES, ANALIZ and CONSIN. 


Figure 2.1 Flow Diagrams For AWALIZ and COPES/COSAIS. 


COPES currently provides four specific capabilities, tuo 
of which will be applied in this work: 
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1. Single analysis - just as if COPES/CONMIN were not 
} used. 

2. Optimization - minimization or saaxiavzization of a aulti- 
wariable function with corresponding constraints. 


COPES requires certain initial data from the user in order 
to coordinate the optimization process. Initial values for 
the design variables as well optimizer control parameters 
are utilized by CONMIN to conduct its numerical analysis. 
There are a faw optimizer parameters which are particularly 
important to the treatment of condenser designs. Ore is the 
finite difference step used in gradient calculations. 
Another is the normalization factor used ir COPES evaluation 
of a constraint function. COPES utilizes the following 
expressions in determining corstraint function violations: 


—RegieBe- 9, 


where SCAL1 and SCAL2 are the normalization factors, BU and 

BL are the upper and lower linits of the constraint, and CFV 

i is the ccnstraint function value. It is intuitive that the 

| normalization factor can play an important role in deter- 

| gaining the size of the active region about a given design 

} constraint. Both finite difference and constraint rnorsali- 
zation will be recalled later during the code development. 

The power of COPES is that it has sinuplified the proce- 

dures involved in using a sophisticated program such as 

COMMIN. The user is therefore fread from the unwanted role 

of systems analyst and can concentrate on the design 

analysis. 


III. CONDENSER DESIGN INPROWESEET PROGRAS (CONDIP) 


A. BACKGROUND 


In the late 1960's, engineers at the Oak Ridge National 
Laboratory developed a sophisticated computer code under 
contract to the Office of Saline Water. This code, called 
ORCON1, [Ref. 2] was generated to aid in the analysis and 
parametric study of large, generally circular condansers. 
Much of ORCON1 was dependent oa Eissenberg's research work 
(Ref. 3] on the effects of condeasate rain on the shell-side 
convective heat transfer coefficient. Johnson [Ref. 5] took 
ORCON1 and made a few minor modifications to determine tube 
side pressure losses and volumetric calculations. Nunn and 
Marto (Ref. 14] further incorporated the correlations 
proposed by Pujii [{Ref. 4} +to determine the effects of 
shearing forces exerted by high vapor velocities on the 
condensate film and resulting shell-side heat transfer 
coefficient. 

It was at this point that the ievelopment of CONDIP was 
bagun. CONDIP was dependent prisarily on the principles 
detailed in ORCON1 but also incorporated subsequent develop- 
ments to the basic prcgram. CONDIP was written , however, in 
such a way as to be ccapatible with tha optimizer, CONMIN. 

CONDIP analyzes a single or double pass, circular or 
senicircular condenser, with steams flowing radially inward 
on the shelleside of the tubes and variable salinity water 
flowing on the tube-side. An optional, rectangular air 
cooler bundle is provwided-for as well as shell-side baffles. 
The circular bundle is normally divided into 30-degree 
sectors with symmetry about the central axis to reduce 
comsputational effort. Unless otherwise specified, tubes are 


21 


placed on a 60-degree eaquilateral, triangular pattern of 
concentric rows with the rows addai from an inaer void ous 
to the outermos* row. The void servas as a collecticn header 
for nen-condensible gases prior to passage througn an air 
cooler, if specified. As in ORCON1, CONDIP proceeds sector 
by sector, row by rew through ths condenser utril’zing an 
average tube to represent the row segment, and calcuiates 
tha following quantities in each sector: 


a) Steam pressure lcsses at the entrance of 2 sector. 

b) Total pressure of the steam/non-condensible gas mixture 
entering a row segment. 

c) Saturation pressure of steam entering a row segment. 

a) Saturation temperature of the steam entering a row seg- 
ment. 

9) Steam flow entering the row segment. 

f) Veloci~y of the steam/non-condensible gas mixture at 
the minimum cross-section in the row segment. 

g) The fraction of non-condsnsible gas in the mixture by 
waight. 

h) Tha overall heat transfer coefficient for the average 
tube in the row segaent. 

i) The steaa-side condansing coefficient. 

j) The tube-side heat transfar coafficient. 

k) The shell-side film heat transfer coefficiant composed 
of the non-condenible gas fila and the condensate fila. 

1) The shell-side friction factor. 

gs) The shellreside pressure loss as steam passes over the 
row sequent. 

n) The shell-side Reynolds ausber based cn the sass flow 
at the sinigua cross-sectional area in the row segment. 

0) The heat transfer rate per square foot 9f condenser 
tube. 

Pp) The sass flow rate of staam/non-condensable gas aixture 
at the sinmua cross-section in the row segment. 
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q) The sass flow of condensate produced as steam passes 
ever the row segment. 

r) Tha cooling water temperature at the outlet end of the 
ccn denser. 

s) The coolant pressure loss on the tube-side. 

+) The average Reynolds number sf the coolant through the 
tube. 

u) The heat transfer coefficient for the non-condensible 
gas fila. 

v) The internal heat film heat transfer coefficient. 

w) The number of tubes per row segment. 

x) The cross-secticn area available for steam flow pez row 
segment. 

y) The cumulative shell-side pressure drop. 

z) The LUTD based cn inlet and outlet coolant temperatures 
and saturation temperature at each row segnext. 

In addition to the above paramatars, the area-weighted 
overall heat transfer coefficients for the condenser, cooler 
and ccmbined condenser are used to calculate the "back 
calculated" log mean temperature difference (LMTD). Steam 
exit-fraction, condenser volume, coolant pumping power and 
numercus other factcrs are calculated from the cumulative 
tasults of the row and sector analysis. 

There are two significant contributions to the external 
fila heat transfer coefficient which have a profound impact 
on the overall analysis. As mentioned earlier, Eissenberg 
(Ref. 3] corrected for condensate inundation effects on the 
external heat transfer coefficient with a series of emper- 
deal relations. He created a flooiing factor F using the 
following relation: 
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where F, is a constant indicating the effect of tube spacing 
and orientation on ccndensate side drainage. With clicsely 
packed tubes, significant side drainage can occur in low 
welocity steaa flow. Condensate generated on tubes above 
may, due to surface tension effects, proceed laterally to 
adjacent tubes rather than down. Thus Py tends to approach 
1.0 for closely packed, staggered tube bundles and zere for 

isperse bundle layouts. S is a constant ctanging in value 


fron: 


{.07 < S <.25) 


If the condensate rain is acting under the influence of 
gravity alone S approaches 0.25. But the influence of any 
steam velocity present begins to alter the rate and dizrec- 
tion ef condensate flow and correspondingly decrease S. Thus 
S is a function of vapor velocity and directon, as well as 
bundle geometry. 

The condensate fila coefficient for the average tube in 
the n-th vertical rcew is then calculated from the uncor- 
rected heat transfer coefficient, hoe as follows: 


J 4 (eqn 3.2) 


i (nF, - (n-1)? - 


n-1 
It is obvious that the determination of a corrected heat 
transfer coefficient is highly dependent on the choice of S 
an Py in equation 3.1. S and Fy are extremely subjective 
constants and there does not exist a current analytical 
expression to detersine then. Yat the choice of these 
constants can have profound impact on the condenser design. 
In CONDIP, as in ORCCN1 [Ref. 2}, the following, relatively 
conservative values fer S and Fy were used: 


S$ = 0.2 


Fa = 0.5 


An additional important contribution tothe ex+ernal 
film coefficient is the effect of velocity shear ferces on 
the condensate film. Pujii [Ref. 4] developed the fcllowing 
experimental correlations to correct the Nusselt number for 
the effect of velocity shear: 

(4a-1) ._ (.5-2a) 

(Nu sya) = ea ae, (eqn 3.3) 
where Nu, is the maan Nusselt number, Rey is the two-phase 
Reynolds number (based on vapor velocity, tube outside dran- 
eter and kinematic viscosity of tha condensate) and Nu, is 


the standard Nusselt number for the zero shear case. The 


empirical constants cand a lie within ths following ranges: 


1.13 < ¢,< 1.24 
0.196 <a < 0.2 


— 


depending on how tube thermal conditions are described. [nr 


CONDIP the values for ¢ and a were: 


c= 1.24 
a = 0.2 


It should be noted that equation 3.3 is only valid in the 


im range: 


| 3.3 < (Re sNu_) < .28. 
L ) 


For sgwaller values of this parameter Fujii recommends the 
ae use of a slightly reduced value of the standard Nusselt 


nuaber: 


Nu = 0.96 wa. (eqn 3.4) 


It is apparent that the vapor velocities conmonly 
encountered in naval condensers can have an impact on the 
heat transfer coeffcient. 


B. CODE DEVELOPMENT 


Johnson [Bef. 5] attempted +9 couple OPCODE2 (his 
version of ORCON1) with the cptimizer CONMIN, but with 
little success. There were several reasons for this. 

ORCON1 uses iterative techniques to solve for such quan- 
tities as condensate rate, steam aass flow rate and steaa 
fressure through the sector. If unrealistic values are 
encountered, such as negative pressures or steam flow, or if 
the final steam exit-fraction exceais a predeterained value, 
ORCON1 stops the analysis, returns to the beginning of the 
program and changes certain initial input parameters. The 
analysis begins again and the process is repeated untila 
satisfactory design is achieved. Thus ORCON1 has a limited 
capability to make design decisions to obtain a feasible 
design. 

CONMIN, as do most optiaizers, requires complete control 
in determining all iterative design variable values. As 
explained earlier, it uses perturbation techniques to calcu- 
lata gradient information for each design variable and 
active design constraint, which it then uses to determine 
search directions. A perturbation of the design varisble by 
CONMIN requires a complete, once-through analysis. %f GRCON't 
is coupled with CONMIN then any adjustment by ORCON1 will 
yield false gradient information to the optimizer and 
hinder, if not completely prevent, CONMIN from arriving at 
the optimum design. During program development it becane 
apparent that the two programs were working independently 
against each other and that in ics prasent state ORCON1 was 
incoapatible with CONMIN. 

In the formulation of CONDIP, it was necessary to locate 
and neutralize all the places where such design decisions 
are nade. By removing the ability for CONDIP to make any 
design decisions it became totally passive and dependent on 
CONMIN for design variable changes. 
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However, once this was accomplished, another problem 


‘area was discovered. In the ORCON! code and subsequently in 


CONDIP there are numerous thermal process and properties 
calculations that use logarithvjic functions and other mathe- 
Batical relationships which could produce singularities if 
the variables in the arguments approach zero or are 
negative. Por example, saturation steam temperature is 
calculated from steam prassure using a logarithmic relation- 
ship. If, during a design analysis, saturation pressure 
approaches a negative value, this raprasents a clear viola- 
tion of physical realities and of the limits of that prop- 
erty. Yet a computer cannot make that distinction so it 
tries to calculate the corresponding saturation temperature 
which, tecause of the logarithaic relationship, would be 
undefined. As just explained, ORCON1 with its built-in deci- 
sion capability simply starts over when this situation is 
encountered. But CCNDIP, being completely dependent on 
CONSIN for design decisions, does not have that capability. 
Remembering that CONMIN requires a complete once-through 
analysis in order to collect enough information to make a 
design decision, it was necessary to somehow bypass such 
mathematical instabilities in ordar to keep the program 
operating. Yet the analysis still had to yield reasonable 
results from the given design in order to obtain meaningful 
gradients. This prompted the formulation of mathematical 
relationships to create “penalty” constraints which, if 
properly written, would indicate to the optimizer that a 
function or thermal property has violated its physical 
limits. However, not only would penalty constraints have to 
te defined, but a "fix up® or "correction" of the violated 
property or variable would be required in order to allow the 
analysis to continue. A good physical understanding of the 
inter-relationship between condenser parameters and the 
thersal precesses resulting from the condenser design is 
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necessary so that the "fix up" of the violated preperty 
would still yield fairly accurate condenser information on 
which CONMIN could base its search for the optimized design. 

For example, a condenser is usually designed around a 
given steam load. If the condenser has too many tubes, is 
too long, or coolant flow is too great, then the condenser 
will be overdesigned. There will be dry tubes within «he 
condenser as all the steam is coniensed before steam flow 
reaches the inner void. In CONDIP this means that zero or 
negative steam flow will be encountered in the analysis. 

If, on the other hand, tube surface area is too small, 
coolant flow is inadequate or the condenser tube spacing is 
too tight and is choking the steas flow, then one of two 
things will happen. Either a quantity of urcondensed stean 
will make it completely through the condenser, or steam 
pressure loss in the condenser will cause steam pressure to 
drop below zero. In addition, there are two reasons why all 
the steam might not condense. It could be simply due to 
insufficient heat transfer surfaca3 area or it could be 
because the saturation temperature 9f ‘the steam has dropped 
below the coolant inlet temperature. If the latter situation 
occurs, then there is no driving force for heat to transfer 
from the steaa to the coolant. Thare is only one way that 
this situation can occur: if steam saturation pressure drops 
bkelow some value indirectly determined by the coolant inlet 
temperature. In aay event this condenser is certainly 
underdesigned and nct capable of supporting the required 
steas load. 

As stated earlier, the purpose of optiaization is to 
obtain the best, feasible design. Thus, an understanding of 
the relationship between the physical characteristics of a 
given design andthe subsequent thermal performances will 
certainly help in defining the appropriate penalty 
constraints and their corresponding liaics. It will also 
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aid in the determination of appropriate “corrections” when 
those limits are violated. It is important to note that with 
the introduction of these penalty constraints, the defini- 
tion of a feasible design is revised. A feasible design is 
now defined as one in which theraal properties and functicns 
are not allowed to violate their physical limits, as well as 
other design constraints, anywhere in the condenser. 

In CONDIP there are three basic thermal properties which 
could create the above mentioned problems if shey fall below 


a certain value. They are steam saturation pressure, steam 
flow and steam temperature. Because of the direct relation- 
ship between steam saturation pressure and temperature, it 
was possible to deal with ther siaultaneously. The solu- 
tions that were developed to overcome the effects of «hese 
thermal violations determined the axtent which CONDIP would 
optinizea. 


1. Steam Flow Effects 


One source of mathematical instabilities within «he 


program is if steam flow over the tubes falls to zero or 
eS below. It is intuitively obvious that steam flow can not 
; physically fall below zero and that in order to keep the 
ie ee program running the steam flow rate aust be kept greater 
te than zero. However, correcting for this alone would 
| : certainly alter the results of that particular condenser 
design, perhaps even imply a feasible design. 

To indicate to the optimizer an infeasible design 
was actually encountered - one in which steam flow had 
dropped below zero - the penalty function WTST was created. 
Since the condenser analysis is performed sector by sector 
and assuming there are J sectors inthe condenser, then 
there had to be arrangements for J penalty constraints. This 
prompted the creation of the array, WTST(J), representing 
constraint penalty functions for eash sector. The absolute 
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magnitude of these functions were iirectly related to both 
the severity 9f the steaa flow violation and the number of 
dry tubes remaining in the sector. WTST(J) ranged in value 
from negative infinity to zero, where a value of zero repre- 
sented a condenser design in which ano flow violations 
occurred. Thus, W8TST(J) were constrained functions whose 
lower limit was zero. For exasple, during an analysis, if 
steam flow was detersined to fall balow zero, then W@TST for 
that sector would be given some negative value. In subse- 
quent designs, as the nuaber of dry ‘tubes approached zero 
and better designs were obtained, then the magnitude of the 
penalty function apfroached zero iniicating no constraint 
wiolation. 

Since penalty constraints are entirely contrived 
trelationships with no real physical basis, it is desirable 
to miniaize their nusber to avoid the possibility of sending 
innacurate signals to the optinizer. 

To eliminate the need to use the WIST penalty func- 
tions as constraints, the values 3f WIST(J) were consoli- 
dated at the end of the sector analysis into the condenser 
steam exit-fraction constraint. Normally, steam exit- 
fraction ranges from zero to one and is simply equal to: 
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By incorporating the WTST violations into steams exit- 
fraction, the exit-fraction was made a continusus function 
Tanging in value from negative infinity ‘to one. The 
negative steam exit-fraction represented a partially dry 
condenser and its sagnitude was in direct proportion to the 
nuaber of dry tubes. Thus instead of having to evaluate and 
calculate gradients for J number of WIST(J) contraints, the 
optiaizer simply had to evaluate 1. previously defined and 
Now expanded exit-fraction constraint. 
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It should be stressed that making steam exit- 
fraction continuous through zero was equally as iaporsant as 
elininating the need for additional constraints. It can be 
reasonably assured that for all practical condenser applica- 
tions, exit steam fraction will always be imited to some 
positive number near zero. flere is where one applies the 
physical knowledge of the condenser and its relation to the 
thermal property of steaa flow. As explained earlier, dry 
tubes represent an ovwer-designed condenser. Thus the natural 
tendency is for the optisizer to alcer those design vari-~ 
ables so as to create a more compact condenser. As this 
occurs, steam exit-fraction will naturally increase. The 
upper active limit of that constraint will determine the 
optinus feasible design. While it is not necessary to have a 
lower limit for steam exit-fraction, it is very important 
for it to be a continuous smooth function especially in the 
tegion near the upper liait. It is therefore critical to 
properly define the fenalty functions WTST(J) in a way so as 
to provide a smooth transition from the negative, artificial 
walues of negative steam exit-fraction to the real, positive 
values. 

Since the steam flow penalty functions will not be 
used as constraints, the analytical results will provide 
gradient information to the optimizer. However, once stean 
flow has been determined to fall below zero, steam flow for 
that first dry row of tubes and all subsequent rows aust be 
fixed up with duamy values to allow the pregram analysis to 
continue. How that "fix-up" is accoaplished will ultimately 
determine the search direction for ths optimizer. 

Physically, once steam flow has gone to zero, there 
should be no further latent heat transfer, no subsequent 
condensate production, and further pressure losses should be 
only due to the flow cf non-condensible gases. It is neces- 
gary to sake the ccaputer generated analysis reflect as 
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Closely as possible these physical realities. Sirce tha 
optimizer no longer has the penalty functions to use in 
calculating a search direction, other constraint values 
obtained from the analysis will dictate the next search 
direction. Gradients will also be calculated using «hese 
results and the determination of the next search direction 
will incorporate these gradients as well. In the case of 
negative steam flow, steam flow and condensate production 
over dry rows were given nominal values which were as small 
as the computer analysis would tolerate. These extresely 
small values closely approximate zero steam flew and 
generate results which reseable physical reality as clesely 
as possible. 

The following example is providad to better illus- 
trate the logic used in CONDIP t> handie negative stean 
flow. CONDIP deterrgines flow rate through each row in each 
sector. During a sector analysis, CONDIP calculates the 
condensate generated at a given row and subtracts that value 
from the steasa flow entering that row to calculate the stean 
leaving. The exiting steams flow rate is then checked to 
detormine whether steam flow has gone to zero. If it has 
not, then the analysis continues. [If it has, then the 
following two events cccur. 

The penalty function, W&TISI(Jj), is calculated for 
that sector and duasy values are inserted for steam flow and 
condensate rate at the row where the violation occurred. 
Por the remainder of the analysis condensate generation and 
subsequent steam flow calculations are bypassed and the 
remaining rows in the sector are fixad up with duany values 
for steaa flow and condensate. The analysis continues 
utilizing these dusay values in all appropriate heat 
transfer and pressure calculations. At the conclusion of 
the sector analysis, the walues of the penalty functions, 
WTST(J), of each sector are incorporated into the stearz 
exit-fraction. 
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If the analysis revealed zero dry tubes ther WTST (J) 
for all sectors would be zero and the stean exit-fracczion 
would sisply be calculated as: 
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If, however, dry tubes were encountered in the condenser 
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analysis then WTS? (J) of some or all the sectors would be 
negative and dependent in sagnituie on tha number of dry 
tubes in each of the J sectors, as well as che severity of 
the steams flow violation. Steam flow leaving any sector 
which has gone dry would be zero and steam exit-fraction 
would be evaluated as: 
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plus a weighted value of all the WISIT(J) penalty function 
values. Using the relationships just described, it is 
apparent that steam exit-fraction: is negative if condenser 
tubes are dary; approaches zero as the design becomes 
feasitle; and is greater than zero if there is steam leaving 
the condenser. 


2. Steams Pressure and Temperature Effects 


. The other possible source of mathematical insta- 
bility occurs when steam pressure falls below some preset 
limit. If pressure falls to zero, humerous aathematical 
singularities will be generated. Yet, before this situation 
can occur steam temperature will have already fallen below 
inlet coolant temperature causing singularities in the log 
mean temperature difference (LMTD) heat transfer calcula- 
tion. Thus, the lower pressure lisit which cannot be physi- 
cally exceeded is not zero but the aininuw saturation 
pressure established by the inlet soolant temperature. In 
CONDIP, this lower limit is given the variable name, PTLIA. 
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As steam flows through the cordenser, pressur3 
continually decreases due to friction losses and therefore, 
it is evaluated at each row in each sector. When the steaa 
saturation pressure drops below PTLIM, indicating a physical 
wiolation of realistic liaits, then the creation of a 
penalty function and a corresponding “fix up" of saturacion 
pressure is required to allow the program to continue. The 
treatment of the problem was therefore analogous to the 
previous situation dealing with negative steams flcw. 

PTST(J) was the penalty function devised +o indicate 
to the cptiaizer that the pressure linit, PTLIA, was 
wiolated in any of the J sectors. Values of thesa 
constraints ranged frem negative infinity to zero, depending 
on the degree and location in the condenser of the viola- 
tion. Since pressure is calculated on a row by row basis in 
each sector, the magnitudes of the pressure penalty func- 
tions were directly dependent not only on how much the 
calculated pressure dropped below PTLIN, but also on the 
number cf rews remaining in the sector. Thus, as the 
condenser approached a feasible desiga thea PTST(J) consraint 
values approached zero, indicating lessening violation of 
the ainiagua pressure. 

As emphasized earlier, it is important to sinimize 
the number of constraints, not only to avoid the possibility 
of sending confusing signals to the optimizer but also to 
reduce cost and iaprove program efficiency. This was accca- 
plished here by inserting dummy values not only into the 
violated pressure variables but also associated thermal 
properties such as condensate generation, heat transfer 
coefficients and heat transfer rates for the row where the 
wiolation occurred and all subsequent rows in the sector. 
The dummy values vere chosen such that realistic gradient 
inforsgation would be sent to the optisizer. The proper 
choice of "fix-up" values for these vwaraiables resulted in 
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the elimination of penalty functions as design constraixts, 
and provided sufficient information to determine subsequent 
search directions. 

It is necessary to understand ‘the influence that 
steam pressure and temperature exert on the overall 
condenser analysis. With this knowladge it will be gasier to 
predict the physical designs which could cause violations of 
the pressure limit. PTLIM is violated due ¢o excessive stean 
flow pressure losses. As explain3sd 3arlier, these large 
pressure losses would result from large steam velocities 
that are found in condensers which ar2 too tightly designed. 
Thus, the particular condenser design is incapable of 
handling the required steam load, implying an infeasible 
design. Understanding this relationship will aid in 
choosing the appropriate "fix-up" values which will indicate 
to the optimizer that when the pressure limit is violated an 
infeasitle condenser has been designed. 

Physically, When steam temperature falls below 
coolant inlet temperature (PTLIM is violated) there is no 
heat transfer from the steam to the coolant and no addi- 
tional steam is condensed. These physical realities must be 
reflected in the condenser analysis. Therefore, in subse~ 
quent rows, condensation and heat transfer rates were set 
equal to zero. Since there is no further condensation, the 
steam exit-fraction is equal to the steaa flow at the point 
of violation divided by the total flow into the condenser. 
Thus PILI indirectly determines the exit steam fraction of 
the infeasible design. This relationship between exit- 
fraction and the PTLIM violation is what makes the penalty 
constraints obsolete. If PTLIMN is violated early in the 
steam's passage through the condenser, steam exit-fractions 
will be large, violating its upper sonstraint liait and thus 
reflecting an underdesigned condenser. As the condenser 
design iaproves, then exit-fractions will decrease. 
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Physically, this can only be accomplished if the condenser 
design "opens up", reducing pressure losses in the 
condenser. Consequently, as conienser designs becom? 
larger, steam exit-fractions decreas2 and the condenser is 
driven tcwards a feasible design. 

The following example illustrates the lcgic employed 
by CONDIP to handle steam pressur2 and temperature viola- 
tions within the analysis. Condenser inlet flow is divided 
by the number of sectcrs in the condenser. Condenser inlet 
saturation pressure is determined by the steam inlet temper- 
ature. Entrances pressure losses are calculated and 
subtracted from the inlet pressure. The resulting pressure 
is checked against PTLIM and a violation at this point indi- 
cates a totally infeasible condenser in which no steam is 
condensed. Steam exit-fraction will thus be equ2l to one. If 
the saturation prassure is graater than PTLIM the analysis 
continues row by row through the sector. Pressure losses 
over each row are calculated and subtracted from the row 
inlet pressure to determine pressure into the next row of 
tubes. If this next-row steam pressur2 is determined to fall 
below PTLIM, then a thermal violation has occurred requiring 
"fix up". Subsequent rows are mads to indicate zero conden- 
sate generation and zero heat transfer. S*team flow over the 
remaining rows is maintained at a constant value, which will 
subsequently be used to determine steam exit-fraction. 
Pressure variables over the regaining rows are given small 
positive values just large enough to allow the analysis to 
continue. Althcugh all heat transfer and condensate caicula- 
tions will be bypassed, the analysis must be allowed to 
continue so that pressure losses will continue to be calcu- 
lated based on the steam flow at the point of violation. 
This is important since stearn flow adjustments to the 
sectors are based on certain pressure comparisons between 
the sectors. The cuargulative sum of all row pressure losses 
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in each of the sectors aust be equal to within some tcler- 
ance. If they are not then steam flow ints each of the 
sectors is altered so that the exit pressures from e?ach 
sectcr converge to scsge coamon value. Thus, anh accurate 
reflection of true pressure losses is important to this 
calculation. 

The value of the steam exit-fraction is again deter- 
mined to be the single constraint necessary to drive subse- 
quent condenser designs to a feasible optimum configuration. 
The pressure penalty constraints proved to be superflous 
information, but the corresponding variable "fix-up" was 
critical in the determination of search direction. 


C. LISITATIONS 


During the develcpment of CONDIP, it became apparent 
that steam exit-fraction would become the key constraint 
during optimization of any objective function. A feasible 
design implies that steam exit-fraction is a small positive 
number perhaps somewhere between zero and 0.1 percent. As 
explained earlier, violaticns of either steam flow or pres- 
sura physical limits resulted in panalty functions and vwari- 
able "fix-up" which were later directly or indirectly 
incorporated in the calculation of steam exit fraction. Thus 
ary feasible design, let alone the optimum one, centers on 
the limits placed upon this design sonstraint. Any number of 
design variable combinations will yield a feasible design, 
and each design variable affects steam exit-fraction differ- 
ently. The intertwined, complex calculations used to ulti- 
mately determining exit-fraction are done by sector and row 
with each design variable repeatedly playing a factor. For 
exanple, the profcund effect of both vapor shear and 
condensate inundation on the shell-side heat transfer 
coefficient and consequently stean exit-fraction, is 
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indirectly determined by numerous dasign variables. However, 
their effects are ispossible to predict. The cascading 
effect of the thousands of calculations performed during +he 
course of a design analysis is to ultimately create a 
Single, highly non-linear variable in the form of the stean 
exit~fraction, upon which design desisions will be made. 

As nore design variabies ware invoived in the analysis, 
the cptimizer had difficulty determining their often 
conflicting effects cn both the objective function and the 
steam exit-fraction. A small perturbation of aach of the 
design variables independently would yield gradients indi- 
cating design improvement. But when these gradients were 
evaluated simultaneously to actually determine the direction 
of the subsequent design, their combined effect would actu- 
ally indicate either no improvement of the objective func- 
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tion or a violation of the steam exit fraction design 
x constraint. The end result would be that the osptimization 
process would stall as no feasible search directicn could be 

obtained. Larger perturbatior’s to the design variables were 

Tequired te properly evaluate their relative effects on the 

objective function and any active or viclated constraints. 

This would enable the optimizer ¢> overcome either small 

inconsistencies or discontinuities in the objective function 

ani the constraint functions which would otherwise prevent 

; the optimizer from reaching the optimum design. This was 
accomplished during data input by changing the normalized 
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finite difference step from 0.01 to about 0.1. Increasing 
the finite difference is not without its drawbacks. As the 
optimum objective value is approached, the optimizer over- 
looks *+he subtle effects of small changes in the design 
warizables because of the relativaly large perturbations. 
Thus, depending on the initial design variables, the 
optimizer will improve the design to some point near, but 
mS necessarily at, the optimun. 
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when a design becomes feasible, st2am exit-fraction will 
always become an active constraint. But the stated goal is 
not in achieving a feasible design but in driving the design 
to a feasible optimus. However, this iterative procass can 
not be accomplished at the expense of violating a constraint 
and it was here that further complication was introduced. 
The initial impetus in any optimization process is to first 
obtain a feasible design. However, onc2 the very small stean 
exit-fractions are obtained «hat are necessary for a 
feasible design, the exit-fraction becomes axtremely sensi- 
tive to any further design variable changes. Thus any effecre 
to further improve the current design could easily cause 
exit fraction to increase. Evan slight increases would be 
perceived as violaticns of the constraint limit and thus 
prevent furthar optimization from the first feasible design. 
There are <¢wo possible solutions to this problem. Bither 
increase the upper limit on the exit-fraction constraint or 
redafine tke constraint. COPES formulates the general 
constraint function in such a way as to allow the user +o 
increase the active region about the constraint limit. This 
is accomplished here by increasing the normalization factor 
in the following expression for tha exit-fraction constraint 
function: 


BU__-~_ _EXITFR 
El seg < 0. 


where BU is the upper constraint bound, EXITFR is the exit- 
fraction constraint value and SCAL1 is the normalization 
factor for this constraint. Increasing the normalization 


factor reduces the optimizer's sansitivity to constraint 
wiolations by enlarging the range of constraint values in 
which the constraint is active. This enlarges the region of 
feasibility and allows the optimizer more flexibilty in 
altering design wariables by reducing the risk of violating 
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the constraint. The cverall effect is that the design opti- 
mization can continue but at the expense of accurate 
constraint limits. The normalization factor used effectively 
for the exit-fraction in CONDIP analysis was approximately 
0.1. 

One of the stated objectives was to create a robust 
program which would consistently yield the single best 
optimum design independent of the initial design and no* get 
hung up on a relative optimum. As it was explained earlier, 
although relative optimums represented design improvement, 
thay also indicated the inability of the optimizer to locate 
the single best or true optimum design. However, the objec- 
tive was achieved for only three design variables. When more 
than three design variables were used, the optimum designs 
becarne loosely dependent on the initial input, although not 
in any predictable way. This is not to say that the 
condenser design did not optimize. By incorporating the 
finite difference and scaling normalization on exit-fraction 
as described above, final designs 1id yield objective func- 
tion values which were continually within about ten percent 
of the true optimum regardless of the initial design. 
However, there was just no guarantee that the single, hest 
optimum design could be consistently obtained. In summary, 
the reasons why CONDIP did not consistently optimize to the 
Single, true optimum were: the extreme non-linearity of the 
steam exit-fraction, the need for a large finite difference 
gradient, andthe need fora normalization factor for the 
exit-fraction upper limit constraint. 

While the optimus design solutions obtained from CONDIP 
may be sufficient, there are several ways to improve the 
results and increase the chances of obtaining the best 
possitle design. The easiest way is to try several initial 
input values until the user is satisfied that the best solu- 
tion has been obtained. The problam with this approach is 
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that it is both costly and time consuming. A second recon- 
mendation is to couple an extremely simplified versicn of 
the condenser analysis with the optimizer t> obtain a 
educated guess as to what the optinum design should be. The 
results of this analysis could thsn be used as inpux for 
CONDIF. OPCODE1, which utilizes the HEI sethods in its anal- 
ysis, is a likely candidate. The aivantage of this approach 
is that a quicker, cheaper analysis can be used to obtain a 
rough idea of the anticipated optimum d2sign. CONDIP can 
then use these design results to obtain aven better and nore 
accurate solutions, faster. There is still no guarantee, 
however, that the true optimum will be solved. Perhaps with 
the development of more robust and versatile cptimsizers, 
ones which uses numerical techniques and methois that ara 
better suited to this type of problem thar CONMIN, sore 
precise solutions can be obtained. However, there is little 
mora that can be done to simplify the analysis of the steaa 
exit-fraction and subsequently linearize the problem. 
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IV. DESCRIPTION OF THE HAIN AND SUPPORTING SUBROUTINES 


A. SAJOB SUBROUTINES 


The following section contains a brief description of 
the sajor subroutines in CONDIP. The appropriate flow 
diagrams are also provided to bettar illustrate and comple- 
ment the explanations. For further information concerning 
the various subroutines and functions see the CONDIP listing 
in Appendix C and ORCON1 (Ref. 2]. 


1. ANALIZ 

This subroutine basically arranges CONDIP in a stan- 
dardized form which is compatible with COPES/CONMIN. CO PES 
uses a variable flag, ICALC, ‘*¢0 coordinate the »ptisization 
process with ANALIZ. Utilizing this flag ANALIZ then calls 
the input, analysis and output portions of CONDIP as 
required. When COPES sets ICALC a3qual to one ANALIZ reads 
in all initial input. This is the only time any input can be 
entered. When ICALC equals two, COPES works with CONMIN to 
optimize the design. ANALIZ makes available the analysis 
portion of the program to be used repeatedly by CONMIN. 
When COPES sets ICALC equal to three, the optimization is 
complete and ANALIZ calls all applicable output subroutines. 
Figure 4.1 illustrates the flow process for ANALIZ. 


2. INPUT 


This subroutine enters alli initial input of data by 
which the initial design is detarsined. The resultant 
design may be either feasible or infeasible, subject to the 
limitations previously discussed, s> it is not critical what 
walues are initially assigned to the design variables. 
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However, the initial input is screanad to pravent the intro- 
duction of totally unrealistic values of variables into the 
program. For exasple, initial tube thickness, tube inner 
diameter, tube number and tube length are all checked to 
ensure that their values are greater than zero. If any of 
the screened initial inputs do not satisfy the aininzua 
requirements, then ¢he program exits prior to entry into the 
optimizer. The limits of the design variales and constraints 
will prevent similar situations from occurring during the 
analysis. Figure 4.2 presents the flow diagram for the INPUT 
subroutine. 


3. ggT3 


This subroutine simply prints all the initial values 
entered in the INPUT subroutine. 


4. QRCON 


This subroutine calculates the bundl2 geometry, 
flooding factors and such coclant flow parameters as pres- 
sure loss, flow rate, and pumping power. Thare are two 
options available to determina bundle geometry, each with 
certain advantages and disadvantagas. 


Option: The number of rows is entered as a constan* and 
the tube number is determined based on pitch, tube outer 
diameter, and row spacing. The advantage of this method to 
determine bundle geometry is that it allows the user to 
linearly vary pitch and/or tube inner diameter by row. The 
disadvantage is that tube number is a dependent variable. 
The optiazer is therefore limited in determining the 
optimua design by the specified number of rows. 


Qption2: The nusber of tubes can be used as a design 
variable while the number of rows is determined by tube 
number, row spacing, pitch and tube outer diameter. There 
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is more flexibility in this method of condenser design but 
it is not possible to linearly vary tube oitch and inner 
diameter. The condenser bundle is generated from a speci- 
fied inner void out, and all the appropriate conderser 
geometry is determined for one of the identical 30 degree 
sectors. Owerall bundle volume is then calculatad 
the ratio of tube hole area to tube sheet area. 


as is 


Once the basic condenser been deter- 
mained, she code then froceeds through an algorithm to calcu- 
late baffle location based on an input value specifying the 
nujaber of baffles desired in the sondanser bundle. After 
this has been completed, flooding factors are determined. 
That is, the number of tubes ina vertical row above the 
central tube in each row is calculatad. don? for 


one side of a 


geometry has 


This is 


each of the six sectors on circular bundle. 


Symmetry is assumed for the other side. These flooding 


Ree, 


factors are later used in calculations +o determine the 
effect of condensate? inundation on shell-side heat transfer 


coefficients. 
: Finally, O&CCN calculates coolant mass flow, coclant 
a velocity, header pressure diffarencs and pumping power based 


. on the type of coolant flow input received. The flow chart 
q in Pigure 4.3 is a simple illustration of the logic used in 
ORCON. From ORCON, the subroutine SECALC is called. 


SECALC 


This subroutine deteraines 


5. 


all the parargeters of 


each of the sectors in the condenser by row. The first 
calculation made in SECALC is the determination of the 
cooler geometry, if there is one. Entrance pressure losses 
. into the condenser bundle are calculated for each sector and 
saturation pressure is checked to ensure that it is greater 
that PTLIS. Prom this point, auch of the reaaining 
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subroutine is comprised of two do-loops with one nested 
inside the other. The outer loop cyclas through however many 
sectors are in the condenser aodel. The inner loop cycles 
through the rows in the sectors. Pressure, tesperature, 
mixture velocity, steam flow and condensate flow are calcu- 
lated at each sector row. The subroutine HETTRN is called 
repeatedly to provide the necessary heat transfer informa- 
tion. Pressure and steam flow is chacked continually at each 
row to ensure that neither falls below its predetermined 
lower limits. In the event that @ither situation occurs, «he 
appropriate penalty function and fix-up procedure is imple- 
mented to enable the analysis to continue. As previously 
discussed, these values are chosen to reflect as accurately 
as possible real conditions which would occur when stean 
flow er pressure violate their physical linits. 

Once all the sectors have been analyzed the cunula- 
tive steas-side pressure losses from each sector are 
compared. Steam pressure at tha inner void must be unifora, 
therefore the sector pressure losses are required to be 
equal within some allowable tolerance. If they are not, then 
the distribution of inlet steam flow into each sector is 
altered to force the pressure loss3s to converge to a single 
value. Once steam flow to the sectors has been adjusted, 
the sector and row analysis in SECALC is repeated until the 
pressure losses within each sector approach a cosgaon value. 
After the pressure ccryparison has been satisfied, certain 
overall condenser parameters are calculated such as steam 
exit-fraction, bundle heat load, and steam-side pressure 
drops. 

Pinally, if a cooler is required, the subroutine 
COOLEX is called. Otherwise the condenser analysis is 
complete. The flow diagram for SECALC is presented in Figure 
4.4, 
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6. HETIEN 


This subroutine is called repeatedly in SECALC to 
solwe for all shell and tube~side heat transfer properties 
for each row in each sector of the condenser. [n particular, 
walues for the owerall heat transfer coefficient and log 
mean temperature difference are utilized by SECALC in 
computing condensate production and heat transfer rate at 
@ach row of tubes. 

On entering this subroutine, a series of estimates 
for certain row wariables are calulated. Based on an assumed 
initial value for the overall heat transfer coefficient, the 
exit coolant temperature and corresponding film temperature 
are calculated. Utilizing these temperatures, the LMTD, 
thermal resistances, individual haat transfer coefficiants 
and numerous other heat transfer parameters are then calcu- 
lated. Pinally, another value for the overall heat transfer 
coefficient is determined based on the above-mentioned anal- 
ysis, and this final value is subsequently compared to the 
initial value. If they are not in. agreement, within a speci- 
fied degree of tolerance, then the initial value for tha 
overall heat transfer coefficient is updated and the entire 
process is repeated until the initial and final values 
converge. This iterative process is necessary as teapera- 
ture dependent heat transfer coeffisients, film temperature 
drops, and exit coolant temperatures are all being calcu- 
lated simultaneously. 

Note that it is in HETTRN that tha concepts of vapor 
shear and condensate inundation are incorporated. Heat 
transfer coefficients are correctei for both effects based 
on the calculations presented earlisar. Also note that since 
steam temperature is never allowed to drop below inilet 
coolant temperature in the calling subroutine SECALC, 
resultant LSTD calculations in HETTRN will not yield 
singularities. 


Once all the heat transfar variablee have been 
determined, control is returned to SECALC where the appro- 


priate results are utilized and stored. The appropriate 
flow diagram for HETTRN presented in Pigure 4.5. 
7. GOOLEX 


This subroutine solves for all th2 necessary parane- 
ters required in the cooler analysis. The cooler is assumed 
to be of rectangular cross-secton with the height of the 
cooler not to exceed the difference between the condenser 
inner and outar radii. The values used for tube pitch and 
tube diameters in the cooler are the same as the innermost 
tow of the condenser bundle. 

Steam exits the condenser bundle, collects in the 
inner woid and enters the bottom row of the cooler. The 
steam then proceeds vertically up through the cooler. The 
physical location of the cooler is not a prerequisite +o the 
subsequent design, although it is expected that the ccoler 
will be placed within the condenser bundle, thus the limit 
on cooler height. 

The first calculation in COOLEX determines the steam 
velocity at minisum cross-secton in the first row of tubes, 
VLCMAX. VLCMAX is directly proportional to the amount of 
steam and non-condensible gas entering the cooler as well as 
the ccoler geometry. Therefore the constrained limits for 
WLCMAX will play a major factor in the overall condenser 
design. 

Subsequent row analysis is treated identically as in 
SECALC. However, all pertinent heat transfer data are calcu- 
lated directly within COOLEZ, making it independent of 
HETTRN. Steam pressure and steam flow are checked at each 
tow to ensure that the appropriate limits are not violated 
and all thermodynasgic parameters are calculated. At the 
conclusion of COOLEX, cooler performance variables such as 
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heat load, exit-fraction steaa, steam pressure losses and 
; overall heat transfer coefficients are calculated and 
_ contrel is returned to SECALC. The flow diagram for the 
COOLEX subroutire is illustrated ian Pigure 4.6. 


8. guTt2 


This subroutine prints the overall condenser bundle 
results including heat load, steaa axit-fraction, overall 
heat transfer coefficient, overall condenser LMTD and bundle 
3 volume. Normally, OUT2 is callei once after the initial 
design is analyzed and again after the optimum design has 
a reen determined. Pinal dasign variable values such as tube 
2 ; number, coolant flow, tube pitch, tube wall thickness and 
tube inner diameter are also printai. 


9. gUur2c 


This subroutine prints the cooler results as well as 
the combined cooler/condenser results. This subroutine is 
called from the subroutina OUT2 and is called only if a 
cooler is required and subsequently designed. Therefore, 
these results will always be printed in conjunction with 
OUT2 output. 
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This subroutine prints a very detailed output of the 
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condenser and cooler results by row and sector. Nearly all 
the thermodynamic and heat transfer proparties are 
presented, thus providing a rather complete picture of 
conditions everywhere in the condenser. This is extremely 
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| : helpful in determining, for examupls, where additional heat 

i transfer enhancement would be most beneficial, or where 

: ‘ baffles should be best located t>. reduce the effects of 
f condensate inundation. 
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B. SUPPORTING SUBROUTINES 


The following is a brief descriftion of supporting functions 
and subroutines called frequently by the main subroutines. 


1. DPSVTY: This subroutine returns the value of the mutual 
diffusivity of the steam and non-condensible gas presert. 


2. XTR: This funcion subroutine transforms the calculated 
data, received in the argument list, to the log values and 
performs a linear regression on two or sore points using 
the aodel. 


3. UPN: This function subroutine calculates the viscosi-~ 
ty of the non-condensible gas in lba/(ft-sec). 


H 4%. BMUFN: This function subroutine calculates the viscosi-~ 

| ty of a saline selution in tha cange of 0-24 percent 

“ concentration and em peratures of 40-210 oF in 
lpas(hr-£ft). 


5. CPAFN: This function subroutina returns a value for the 
heat capacity of the inert, neon-condensible gas mixed in 
with the staaa in units of Btu/ (lba-n0l-OF). 


} 6. CPFN: This function subroutina calculates the specific 
1 heat of a saline sclution units of Btu/(lbm-°F). 


7. CPSPN: This fuaction subroutine calculates the heat ca- 


oF 

+ pacity of stearz in Btu/(lba-201-9F). 

| : 8. HPGFN: This function subroutine returns a value for the 
;: latent heat of vaporization of water in Btu/lba. 
: 9. PRSDRP: This subroutine returns the shell-side pressure 
drop across a row of tubes in psia. 
z 


10. PSATPN: This function subroutine calculates saturation 


pressure of steaz as functicn of temperature. Pressure is 
returned in units of psia. 


11. BOEFN: This function subroutina calculates the censity 
of a saline solution ofconcentration range 0-24 percent and 
temperature range of 40-300 OF. Dansity is returned ir 
units of lba/(cu.ft.). 


12. SKBFN: This function subroutine calculates the thermal 
conductivity of a saline solution of concentration range 

0-24 percent and a temperature range of 40-300 °F. Thermal | 
conductivity is in (Btu)/(hr-f-9F). 


13. JSATPN: This function suroutine returns the value for 
st@gam temperature in °R given a pr2ssure in psia. 


14. VGFN: This function subroutine calculates the specific 
volume of steam as a function of temperature and pressure. 
It has units of (cu.ft.)/lbn. 


15. SWITCH: This function subroutine reverses the order of 


a stored array. 
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Figure 4.1 Plow Diagrarm for the ANALIZ Subroutine. 
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Figure 4.2 


Flow Diagram for the INPUT Subroutine. 
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ORCON Subroutine. 


Plow Diagrarg for the 
53 


CALCULATE 
PLOODISG FACTORS 


Figure 4,3 


B CO) LER 


1] 
L 
GEOMETRY 


CALCULATE 
GULA 


RECTAN 


00 Loop 
EctTorRs ti 


FOR S 


er ec © mm napa me ei 


for the SECALC Subroutine. 
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Diagram (continued) 
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Plow Diagrasa for the HETTRE Subroutine. 
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Figure 4.5 
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Flow Diagram for the COOLEX Subroutine. 
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COOLEX Flow Diagraa (continued) 
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A. CONDIP VERIFICATION 


It was desirable to verify the single pass capability of 
CONDIP (i.€. without optimization) as a predictor of 
condenser perforagance by comparison with actual experimental 
data. However, ccmplete and accurate data on condenser 
design and corresponding performance is not always readily 
available. Lynch (Ref. 18] encount2red this same prceblem in 
attempting to verify ORCON1. However, he did marage to 
locate scme actual experimental cordenser data, obtained 
during a test conducted to determina the general performance 
of the DDG-37 class propulsion machinery [(Ref. 16]. The 
test took place at the Naval Boiler and Turbine Laboratory 
and was conducted primarily to determine the performance of 
the turbine and reduction gears. Some limited condenser data 
was taken as a by-product. The various measurements were 
obtained as described below: 


1. Steam flow measurements were made by weighing «he 


condensate. 


2. Cooling water inlet and outlet temperatures were meas- 


ured by thermometers installed in the inlet and discharge 
lines. 


3. The heat load is calculated based on total steam flow 
into the condenser sultiplied by the difference between 
inlet steam and condensate enthalpias. 


4 Circulating water flow was determined from a heat 
balance around the condenser. The total heat load was 
divided by the circulating water heat capacity and 
temperature rise to cbtain flow rate. 
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5. Condenser inlet pressure was determined by pressure 
instruments located above the condensar inlet flanges. 


6. Non-condensible gas flow was maasured by a flowrator. 


7. Pressure at the air ejectcr was measured directly. This 
pressure and the inlet pressure determined the pressure drop 
across the condenser. 


It should be pointed out that this data were not recorded 
with the care that normally accompanies scientific data 
collection. Neither the instrumants nor the techniques 
eaployed were particularly accurate. The possibility that 
this observed data are in error casts a cloud over the cred- 
ibility of the corresponding condensez performance, which 
was calculated based on those values. However, for a lack of 
better alternatives, this data and the resulting condenser 
analysis will be used to determine the reliability of 
CONDIP. 

The DDG-37 condenser geometric design variables obtained 
from the technical manual [Ref. 15] and input parameters 
corresponding +o a full speed run ara presented in Table I. 
An attempt was made to repeat the design using CONDIP. The 
results cf CONDIP*s proposed design as well as the experi- 
mental performance are presented for comparison in Table II. 
Percentages were calculated to quantify the differences 
between the actual and hypothetical performance. 

Before elaborating on the results of ¢his verificaticn, 


some notable differences between the two designs must be 


clarified. First, a specific fouling factor was not deter- 
mined at the time of the experiment and was therefore not 
provided. A somewhat realistic claanliness factor of 87.5 
percent (fouling factor of .0002) was utilized. 

Second, in the DDG-37 condenser the rectangular cooler 
appears to te inserted directly into the condenser bundle. 


However, in order to accomodate the cooler, the bundl2 nus+ 
expand or distort. In addition, a void of some dimensions 
4 must be provided-for in the center of the bundle tc collec 
any uncondensed steam and non-condensible gases. Diagrams in 
{Ref. 15] indicate that the DDG-37 condenser is indeed 
nearly elliptical in shape with bundle axes of 5.67 and 7.17 
feet. Although it is apparent that a void does exist, exact 
@imensions can not be readily determined frem the available 
information. 

CONDIP approximates this design by creating separately 2 
circular condenser bundle anda ractangular cooler, the 
height of which cannot exceed the difference between the 
outer and inner bundle radii. A circular void of pre- 
determined size is provided-fer when determining the 
condenser bundle geosetry. Subsequent volume calculations 
are performed on the condenser bunile and cooler separately 
and the overall condenser volume is computed as simply their 
sum. 

Although CONDIP dces not exactly duplicate the geometric | 
configuration of the DDG-37 condenser, it was possible to 
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Banipulate cartain initial design variables in order to 
cause CONDIP to develop an approximately equivalent configu- 
ration. These variables were chosen because, for saall 
! changes in their values, there is a rather significant 
! change in the bundle geometry with relatively snall effects 
on the overall condenser performance. Since there are no 
4 specific dimensions provided for thea inner void in *he 
. DDG-37 technical manual, it was picked to be one of the 
design variables to be adjusted. Row spacing was also 
adjusted because it satisfied the conditions described 
above. Through trial and error a combination of row spacing 
and inner void radius were determined, froma which CONDIP 
yielded a geometric design similar to the DDG-37 condenser 


and that satisfied condenser requiresaents specified in 
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(Ref. 17). In this particular case, the void diameter ard 
row spacing were determined to be 1.1 feet and 1.35 inches 
respectively. This arrangement enabl2ad the condenser model 
to closely approximate the tube sheet area ratio of the 
actual DDG-37 condenser. This manipulation, however, must 
be interpreted as another source sf error and innacuracy 
when comparing CONDIP's conienser performance with the 
experimental results. 

Lastly, condenser designs in (Ref. 15] reveal that chree 
different «ube patterns were employed in the DDG~-37 
condenser. [In addition, two differant values for tube pitch 
were used - a pitch of 1.4 in the condenser bundle anda 


pitch of 1.3 in the ccoler. This situation cannot be dupli- 
cated in CONDIP. Therefore a constant pitch of 1.40 and a 
uniform tube pattern were utilized throughout the condenser. 

The design approximations utilized in CONDIP to try to 
geometrically simulate the actual DDG-37 condenser intreduce 
significant uncertainty into subsequent design comparisons. 
This, coupled with the fact that the data collected is also 
suspect, would imply that it is rather difficult to verify 
CONDIP's analysis with the information available. + should 
also be noted that CONDIP is sensitive to even small vari- 
anc3 in either the data collected (i.e. steam inlet tempera- 
ture) or the approximated design variables. However, 
despite the above-mentioned problems associated with 
equating CONDIP*s condenser to the DDG-37 condenser, the 
experimental data obtained from tha DDG-37 condeaser still 
provide the best available base upon which to make a reason- 
able determination of CONDIP's capabilities and limitatioas. 

In comparing the results in Table II, it is immediately 
clear that there is significant diffezcence between certain 
condenser performance parameters pradicted by CONDIP and the 
corresponding experimentally derived condenser perfcrmaace. 
Alraady, much has been said about the numerous geometric 
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approximaticns used to model the DDG-37 condenser. Buc 

questionable data and geometric Manipulations may rot 

completely explain the 13 percent axit-fraction and gerezral 

, poor performance generated in CONDIP*s analysis. I+s values 

for the average overall heat transfer coefficient and heat 

rejected were significantly lower than the 2xperimental 

results. One source of the problem lies in the actual heat 

transfer analysis performed in th2 code. Lynch [Ref. 18} 

graphically illustrated how sensitive this analysis is to 

the effects of condensate inundation. In particular, by 

making small changes - within the allowable ranges - in the 

constants used in Eissenberg's correlations for condensate 

rain, Significant improvement could be crealized in «he 

overall heat transfer characteristics of the condenser. 

i CONDIP*s results, when compared t2 the experinarntal data 

support the argument that the values currently used in the 

inundaticn correlations are rather conservative in nature, 

and cause the overall analysis to yield a poor preformance 
for the given steam load and condenser design. 

Therefore, in order to present CONDIP with a fair test 
to determine its credibility as a design pradictor, some 
additicnal work must be first accomplished. A condenser 
geometrically identical to the ga2neral model created in 
CONDIP should be constructed with complete and accurate data 
acquisition systems tc establish a thorough data base from 
which to compare. Also, more research should be performed on 
i the effects of condensate inundation and velocity shear to 
t obtain more precise correlations in determining their 

overall effects on the film heat transfer coefficients. 

One last additional point should be mentioned. In 
comparing the steam-side pressure drops through the 
condenser, it was shcwn that CONDIP'’s pressure losses were 
nearly 72 percent larger the the actual physical measure- 


ments. However this radical differance is mainly due +o the 
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high pressure losses experienced at the entrance of the 
cooler as alarge volume of steam tried to force its way 
through the small available area. Therefore, the signifi- 
cance of this large disagreement in results is relatively 
minor and can be treated siuply as a consequence of the more 
important heat transfer limitations in tha comparison run. 

Although the goal of verifying CONDIP as ae design 
predictor has proven elusive, it was still possible to 
demonstrate its capabilities through comparison studies. 
Therefore, the remaining emphasis in this thesis is to show 
the ability of CONDIP (in combination with the optimizer 
COPES/CONMIN) to take an initial dasign with a given frame- 
work of constraints and design variables, and obtain better 
designs based on a desired objective function. 


B. EXPLANATION OF THE CASE STUDIES 


The following case studies were devised to best exercise 
the capabilities of CONDIP. They were made as realistic as 
possible so as ¢t9 simulate the problem of condenser design 
and specification ccenfronting the engineer during the early 
stages of power plant design. The condenser performance 
returned by CONDIP during the verification run and contained 
in Table II will serve as a baseline for comparing the 
results cf each case study. The baseline condenser perform- 
ance is based on the design parameters from the DDG-37 
condenser listed in Table I. It was stated 2arlier that 
CONDIP*s optimization results are slightly sensitive to the 
initial design if more than three design variables are used. 
Since all the cases involve eight or more design variables 
it would be best, for the purposes of comparison, to start 
from the same initial design in all cases. Therefore, the 
initial design variables used for the verification cun and 
contained in Table I will be utilized as the baseline 
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design. Although many of these initial design variables will 
be allowed to change during optiaization, certain tasic 
condenser requirements will not. They include: steam flow 
into the condenser, inlet steam saturation pressure and 
teaperature, cooling water injection temperature, the frac- 
tion of non-condensible gases in the steam, the tube fouling 
factor, and the tube material. It should be noted that 
although there was an initial value for row spacing given in 
Table I, row spacing was not used as a design variable 
during any of the optimizations. Instead, the program used 
the default method of row spacing calculation available in 
the ccde where the rows are spacei such that a 60-degree 
equilateral triangle pattern of concentric rows is obtained. 
Row spacing is therefore dependent on tube pitch and tube 
outer diameter by the following relation: 


RSPA = (SDDO * ODOT) * .866 (eqn 5.1) 


where &SPA is row spacing, SDDO is tube pitch, and ODOI is 
tube outer diameter in inches. 

There are a few key points t> be kept ia mind when 
comparing the results of the case studies with the baseline. 
Pirst, the baseline design is an infeasible and inadaquate 
design. Its performance indicates that it is not capable of 
Supporting the required steam load by returning a_ stean 
exit-fracticn in excess of 13 percent. So any gains in the 
objective function that were realized in the case studies is 
even more remarkable since it is a necessary condition that 
the optimua design be a feasible dasign, defined as having 
an exit-fraction not greater than 1 percent. Second, the 
percent change referred to when analyzing the results is 
calculated based on the baseline design. Thus the baseline 
serves as a unifora frame of reference. Next, it should be 
noted that because of the large nuaber of design variables 
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and constraints, intuition on how an optimized result will 
turn out is not always applicable. Pinally, it wiil be 
easier to understand the effects of the vwarious design 
parameters by keeping in mind the following, very basic, 
heat transfer correlation: 


Q20 * A * LSTD (eqn 5.2) 


where Q is the rate of heat released as the steam condanses; 
U is the overall heat transfer coefficient; A is the heat 
transfer surface area; and LMTD can be interpreted as the 
thermal driving force between the steam and the coolant. Q 
is directly dependent on steam flow and pressure into the 
condenser, the percentage of that steam that is condensed, 
and any subcooling of the condensate. Por a given steam load 
and a very small exit-fraction Q is nearly constart as the 
optimized results in all the case studies will indicate. 


1. Constraint Framework for CONDIP 


In order to simulate an actual trade-off study, the 
constraints and their respective limits were kept constant 
for all the case studies. The conienser was to be designed 
with a maximum bundle diameter of ten feet, a maximum and 
minimum tube outer diameter of 0.625 and 1.25 inches respec- 
tively, a steam exit-fraction of not more than 1 percent, a 
maximum cooler inlet velocity (VLCNAX) of 200 feet/second, 
and a ratio of tube sheet hole area to total tube sheet area 
of less than 0.30. 

The constraint on bundle diameter was chosen sone- 
what arbitrarily. It seems unlikely that this limit would be 
realistically exceeded, although certainly space require- 
ments would dictate the exact configuration. Tube outer 
diameter is dependent on the values for tube wall thickness 
and tube inner diameter. Thus the liaits imposed on tube 
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outer diameter represent realistic crestrictions on the 
possible cobinations of inner diamater and wall thickness. 
These restrictions are based loosaly on anticipated tube 
structural and strength requiremants and correspond to 
values of normally available tubes (Ref. 19]. The maxinua 
limit of 200 feet/second for VLCMAX was also a somewhat 
arbitrary but realistic limit. It is assumed that stean 
velocities cften exceed that value in the condenser bundle. 
It is recalled that steam axit-fraction will play a 
Fs significant role in the determination of the final optinun 
* design. The baseline exit-fraction of 13 percent predicted 
a ; by CONDIP for the DDG-37 condenser is unsatisfactory. 
4 Tharefore a more reasonable upper limit of 1 percant was 
placed on this constraint. Although CONDIP will return a 
Buch more ccnservative design if 1 percent vice 13 percent 
- is used as the upper limit, the subsequent design will be 
F such more credible. 
“ Finally, the amount of tube sheet material that can 
q be reaoved by drilling for the installation of condenser 
tubes is. specified at 24 percent of the total tube sheet 
area in (Ref. 17). This area ratio imit represents a 
structural liait imposed to ensura that the tube sheets do 
not fail due to heat and pressure stresses in the condenser. 
However, CONDIP does not take into account the space between 
the condenser and tube shell normally used in area ratio 
calculations as blank tube sheet area. For this treason and 


: ‘ to allow mora flexibility .in the design analysis, the 
me) 7 constraint limit was set at 30 percant. 
\ In summary, the general dasign constraints and the 
associated upper and lower bounds ware: 
i 0.625 $ tube outer diameter (inch) $< 1.25 
: ; 1.0 $ bundle diameter (feet) < 10.0 
a steam exit-fraction (%) < 1.0 
i VLCHAX (ft/sec) S$ 200.0 


area ratio < .30 
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7 These design cecnstraints ani associated bounds were 
used inall the case studies except wher? specifically 
modified. | 


4 2. Resign Variable Framework for CONDIP 


At least eight design variables were usad in all the 
case studies. They include tube inner diameter, tube wall 
thickness, tube pitch, the number of tubes in the condenser, 
tube length, the inner void radius, tha percent of the tubes 
in the cooler, and cooling water velocity. Side constrainrts | 

i were placed on all cf these variables t9 correspond to | 

;* : either realistic physical limits or available standardized | 

7 materials. 

: Tube wall thickness was not allowed tos fall below 
0.022 inches (BWG 24) or exceed .109' inches (12 BWG), sizes 
normally available commercially. Tube inner diameter was 


' restricted to values between .407 and 1.206 inches so as to 
aq yield tube outer diameters, within the limits specified 
- earlier. 

‘ Tube pitch is defined as tha ratio of the center to 


center spacing between adjacent tubes in a row to the tube 
outer diameter. Tube pitch is an accurate measure of how 
closely packed the tube bundle is. Generally accepted 
values for pitch lie in the range of 1.3 to 1.7. However, 
to provide more latitude in the design process this design 
wariable was allowed to wary in the range between 1.1 to 
2.0. 
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There is no guidance availabla as to the allowable 
range for tube length in the condenser. Since the lower 
ligsit was not expected to be crucial, it was set randomly at 
1.0 feet. The upper limit of 25.0 faet was a realistic lisit 

§ considering the size of the tube diameter being worked with. 
Inner void radius and the percent of tubes in the cooler 
were chosen to be design variablas simply to enhance the 
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flexibility of the code in designing the conderser model. 

The bounds for both variables wera entirely arbitrary with 

only ccmmon sense as the determining factor. The upper and 

lower limits on the percent of tubes in the coola= was 
_ established as 10.0 and 2.0 percent raspectively. The upper 
by and lower bounds on the inner void radius was set a+ 1.0 and 
pe 0.1 feet. 

Cooling water velocity generally ranges from three 
to nine feet per second in value for all common tube 
materials, except titanium which has an upper limit of 15 
feet per second. Exceeding these upper limits risks exces- 
Sive tube erosion and material damage. Finally the number of 
tubes was perzitted to vary between 1000 and 8000 tubes for 
the purpese of improving design flaxibilty. It is extremely 
unlikely that, for most propulsion applications, tube number 
would fall below 1000. The upper limit was simply chosen as 
j a realistic cutoff point in terms of complexity, cost and 
of maintainability. 
In summary, the general design variables and the 

associated side-constraints were: 


0.407 < tube inner diameter (inches) < 1.206 

0.022 s tube thickness (inches) < 0.109 
| 2.0 $ percent of tubes in cooler < 10.0 
0.10 < inner woid radius (f2et) < 1.0 
3.0 S$ coolant velocity (ft/sec) < 9.0 

1.0 $ tube length (feet) < 25.0 
1000 < tube number <¢ 8000 
1.1 $ tabe pitch < 2.0 


As a, 


As in the case of design constraints, these design 
variables and their respective lisits were used consistently 
in all the case studies unless otherwise specified. 
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C. CASE STUDIES USING CONDIP 


1. Case One 


The objective of this case was to minisaize condenser 
volume. The final results of the optimization along with the 
initial parameters is listed in Table III. 

These results show a 16 percent decrease in 
condenser volume with a corresponding 24 percent increase 
in pumping power. The source cf tha improvemen* can be 
understood by noting the following: 

1) Tube wall thickness was reduced from 0.049 to 0.022, the 
minimum side-constraint, thus allowing tube inner diameter 
to incraase while aaintaining a minimum tube outer diameter. 

2) The number of tubes shrank slightly as did tube length, 
resulting in a smaller heat transfar surface area. 

3) Tube pitch increased aarkedly, causing a reducticr in 
steam pressure losses which then ansured that high values 
for steam saturation pressure and tamperature would be sain~ 
tained throughout the condenser. fhe large pitch aiso 
reduced steam velocities, allowing the cooler inlet velocity 
limit tc be satisfied. Row spasing decreased from «he 
initial value of 1.35 inches, thus decreasing condenser 
volume. 

4) Cooling water velocity increased to the saxisgua allo- 
wable value of 9 ft/sec which correspondingly resulted in 
larger head losses and coolant flow, causing overall pumping 
power to increase. 

As cooling water velocity increased and tube wall 
thickness decreased, then their raspective thermal resis- 
tances were disinished. The cumulative effect was to 
improve the overall heat transfer coefficient. LATD rose 
primarily asa result of the higher steam teuperatures 
throughout the condenser. It is apparent by looking at 
equation 5.2 that increasing the driving forces for heat 


transfer, such as the overall heat transfer coefficient anid 
LMTD, allows the heat transfer surface area to decrease. 
This resulted in a similar reduction in condenser vclume. 

The constraint limits that prevented further design 
improvement were the upper bound on the cooling water 
welocity, the upper limit on the tube sheet area and the 
upper limit on VLCMAX. 


2. Case Ivo 


, The objective of this case was to minimize the 

puaping power required to overcome the «ube-side head losses 

u and drive the cooling water through the condenser tubes. The 

final results of the optimization are presented along with 
the initial design in Table IV. 

The results indicate a dramatic 90 percent reduction 


ee 


in reguired pumping power with an 2qually large 120 percent 
increase in condenser volume. The aajor factors involved in 
the design improvement along with their relative effects are 
briefly explained below: 

1) Tube inner diameter increased 27 percent while tube 
thickness remained relatively unchanged. Thus tube cuter 
diameter was caused to increase. 

2) The number of tubes in the condansar rose significantly, 
along with tube length. This, coupled with the enlarged tube 
outer diameter resulted in nearly doubling the haat transfer 

area. 

‘ 3) Tube pitch icreased 29 percent, which allowad stean 
. saturation pressure and temperature to be maintained at 
consistently large values inthe condenser. This had a 
benefiting effect on the associatei LHTD calculation. The 
large pitch also helped satisfy the steam velocity lisit 
into the cooler. The tube spacing dacraased from the initial 
walue of 1.35, but by a smaller amount than the previous 
case because of the large values for tube pitch and outer 
diasete:. 
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4&) Cooling water velocity dropped to the minimum allowable 
liait of 3 ft/sec. This had tha effact of reducing tube-side 
head losses and ccolant flow through the condenser. 
Consequently, pumping power was drastically reduced. 

The combined effect of all these charges can again 
be put in perspective by looking at equation 5.2. For the 
given steam and corresponding heat load, the heat transfer 
area increased drastically, allowing both LMTD and the 
overall heat transfer coefficiant to decrease. A smaller 
overall heat transfer implies a smaller convective tube-side 
contribution which in turn permits coolant velocity to 
reduce to its lowest allowable valus. The LMTD decrease i 
explained by the fact that cooling water was spending sore 
time in the tubes, thus causing the average cooling water 
temperature to rise. However, tha subsequent reduction in 
LMTD was minimized by the fact that a high steam temperature 
was maintained in the condenser. 

There were no active constraints in this design 
outside of cooling water velocity which prevented further 
design improvement. However, the penalty paid in terms of a 
huge condenser volume, appears prohibitive. 


3- Cass Three 


The objective of this case was to miniaize pumping 
power while holding condenser volume constant a* the initial 
value of 432 cubic feet. This was a particularly interesting 
test case as the results in Table V bear out. The required 
pumping power was reduced by nearly 38 percent with no 
change in volume. The effects of the design changes which 
resulted in the design improveasn* are presented below: 


1) Tube inner diameter increased noticeably. However, the 
effects of this increase on tube outer diameter was sini- 
mized by a large drop in tube wall thickness. Thus, tube 
outer diameter remained relatively unchanged. 
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2) The number of tubes experienced a minor reduction, while 
tube length increased. The overall affect was to increase 
heat transfer surface area. 

3) Tube pitch again rose by nearly 25 percent, causing 
steam Saturation temperature and pressure ¢9 maintain 32 
nearly constant value throughout tha condenser. This hada 
beneficial effect on the LMTD between the steam and the 
cooling water. The larger pitch also had the additional 
effect of reducing steam velocity thus allowing the subse-~ 
quent design to satisfy the upper limit on steam velocity 
into the coclar (VLCMAX). The combination of tube pitch and 
tube outer diametar resulted in a reduction in row spacing 
from the initial value of 1.35 

4) Cooling water velocity decreas2d by about 21 percent. 
This effect was manifested in subsequent pressure head, 
coolant flow and pumping power calculations. 

Looking at equation 5.2 we see the same general 
pattern emerging as in Case 3, but with more subtlety ir the 
changes. Heat transfer increased, but not at the expense of 
volume, Cooling water velocity was allowed to decrease 
while the cverall heat transfer coafficient actually rose. 
One explanation is that as the tube wall got thinner its 
thermal resistance got smaller which more than offset the 
loss of convective heat transfer contribution from the 
coolant. The LMTD dropped slightly iue to the higher average 
coolant temperature of the coolant in the tubes. 

The constraints which becaze active and prevented 
further improvement in the design include tube sheet area 
ratio as well as tube wall thickness. However, tube wall 
thickness was particular crucial because of its related 
effect on heat transfer. 
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4. Case Four 


The objective of this case was to minimize condenser 
volume while holding pumping power constant at the inicial 
walue of 55.7 horsepower. The results of this case can be 
found in Table VI. Chesen *o contrast the results in Case 4, 
tha relative improvement in this design objective was not 
nearly so impressive. Condenser volume shrank by only 12 
percent. An explanation of the causes and effects is 
provided below: 

1) Tube inner diameter increasei, with 2 corresponding 
dacrease in tube wall thickness to yield thea minimum allo- 
wable tube outer diameter. 

2) The number of tubes decreased noticeably, tending mini- 
mize bundle volume. Note, there was only slight increase in 
tube length. The overall effect was to similarly reduce heat 
transfer surface area as condenser volume decreased. 

3) Tube pitch again increased significantly, having the 
same effects on steam pressure, tamperature and steam 
velocity into the cooler as discussed 2arlier. A large tube 
pitch benefits the LMTD between the steam and the cooling 
water. Row spacing was again a factor in reducing condenser 
volume as before. 

4) Cooling velocity decreased slightly as did head loss. 
But overall coolant flow increased due to an increase in 
tube inner diameter. The net effect was to maintain pumping 
power. 

Again, referring to equation 5.2 , it is clear that 
the slight decrease in heat transfar area was offset by the 
Slight rise in LSTD resulting froa higher condenser stean 
temperatures. The significant improvement in overall heat 


transfer coefficient, therefore, is what makes the heat 
balance work. The large decrease in tube wall thickness and 
corresponding reduction in thermal rasistance contributed 
heavily to this imprcvesent. 
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There were several contraint limits which prevented 
any additional objective optimization. They include the 
minimum tube wall thickness, tube sheet area ratio and stean 
velocity entering the cooler. 


The objective of this case was to minimize condenser 
volume while exercising CONDIP's csapabilities to linearly 
vary tube pitch and tube inner diamatar by row. Thirty-five 
rows were used, which was the identical number as ths 
initial design. Tube pitch and tupe inner diameter of both 
the outermost and innermost rows ssrvad as design variables 
in this case. However - because tupe number is now a depen- 
dent variable based on the number of rows, tube pitch, and 
tube diameter ~- it cculd net be used as a design variable. 
The optimized results of this analysis along with the 
initial design are presented in Table VII. 

The results of this test cas2 indicate a condenser 
volume which is 20 percent smaller than the initial design 
as compared to a 16 percent decrease in Case 1. The basic 
reasons and axplanations as to why volume was able to be 
reduced remain fundamentally the same as in Case 1, 
Attention will therefore be focussed on the effects of 
linearly varying pitch and tube dianetar. Final tube pitch 
ranged in value from 1.75 in the inner row to 1.44 in the 
outer row. Similarly, tube inner diameter ranged from .729 
in the inner row to .583 in the outar row. 

It is believed that smaller pitch and inner diameter 
were used in the outer row because of the higher available 
steam saturation pressure and temperature. The resulting 
higher steam velocities enhanced the beneficial effects of 
wapor shear on the external heat transfer coefficient 
thereby improving heat transfer on the outer rows. As stean 
pressure decreased, then tube pitch and tube inner diameter 
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increased to compensate and extract all the available heat 
from the steam. Consequently, steam velocity decreased and 
was able to satisfy to the limit imposed on the ccoler 
entrance velocity. The end result is a condenser geomexsr 

that makes complete use of “he available resources and 
conforms tc the geometry to take advantage of the thermal 
conditions in the condenser. The big limitation with this 
approach is that the number of rows is held constan*. Thus 
the subsequent condenser is designed around that value and 
the subsequent optimum design is a function of the number of 
rows specified. 
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TABLE I 
Input Design Data 


Panne ee Reece 


are a een nents eo 
PARAMETER | VALUE 

A (cas! 
Total number of tubes { 5230 | 
{ Tube length (feet) 19.3 {| 
{ Tube inner diameter (inches) 9.527 | 
| Tube wall thickness (inches) 9. 049 | 
Tube outer diameter (inches) { 0.625 
| Tube mat'l thermal conductivity { 25.0 | 
btu/ (ft-hr-9 F) 
Tube pitch ( 1. 38 { 
| Percent of tubes in the cooler { 7.0 
Steam inlet flow (lbm/hr) 16 1,961 
Fraction of non-condensible gas 37.1 | 
: | 5 | | 
‘ Steam inlet pressure (psia) 1. 294 ‘ 
| Steam inlet temperature (9F) 110.52 | 
} Coolant inlet velocity (ft/sec) | 8.473 
: | Coolant inlet temperature (°F) { 75.66 | 
a Fouling factor -0002 
q ; \ Inner void diameter (feet) { 1.1 

| Row spacing (inches) | 1.35 
[ _ eee eee enrenae apy aeennees emen eens mee — | 
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CONDIP Verification Results { 


ee a ne ne a ny 
| TABLE II 
{ 


EXPERIMENT | CONDIP CHANGE | 
FARAMETER RESOLTS RESULTS (%) 

| { 
Heat transfer area 8805 i 8814 +0.10 

(sq. ft.) | 
Overall heat transfer 635.2 { 547.9 | -9.5 | 
coefficient | | { 
btu/ (hr-sq. ft.-9F) | | | 
Log bean temperature 28.24 1 28.62 |{ #1.3 { 
{ aifferenca (FP) { | { { 
| Coolant temperature | 10.61 H 9.90 -6.7 | 
ae | | ee 
Coolant mass flow 1.503 1. 540 | 42.5 | 
rate (107 lbm/hr) | | 
ioe volume | 2<---- | 432.3 | -o--- 

cu. ft. 

{ { 
| Bundle diameter (ft) 5.7 7.2 7.17 | worn 
Shell-side pressure | 0.751 1. 29 +71.5 

drop (psia) 

Stean gxit- fraction aatateeatal 13.3 -o--- 

(% of input) 

Hea+ rejected 1.595 { 1.451 -9.03 | 
\ (108 btu/hr) | | 
\ Area ratio 0.291 H 0. 266 | 78.59 
= | ees oe 
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TABLE III 
Volurgae Miniaization 


PARAMETER 


Total number of tubes 
% of tubes in cocler 
Tuke length (ft) 

Tube inner diam. (in) 
Tuke wall thick. (in) 
Tube outer diam. (in) 
Tube pitch 

Void diameter (ft) 
Bundle diameter (ft) 


Condenser volune 
(cu. ft.) 


Area ratio 


Coolant inlet vel. 
(ft./sec) 


Coolant mass flow 
rate (107 lbm/hr) 


Head loss (ft H20) 
Pumping power (hp) 


Coolant temperature 
tise (°F) 


Log méan temperature 
ference (°F) 


Heat transfer area 
(sq. ft.) 


Average overall heat 
transfer coefficient 
btu/ (hr-sq.ft.- 


Steam exit-fraction 
of input) 

£ Py a 

® bta/hr) 
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TABLE IV 
Power Minirnization 


PARAMETER 


Total nuaber of tubes 
% of tubes in cocler 
Tube length (ft) 

Tube inner diam. (in) 
Tube wall thick. (in) 
Tube outer diam. (in) 
Tube pitch 

Void diageter (ft) 
Bundle diameter (ft) 


Cendenser volume 
(cu. ft.) 


Area ratio 


Ceolant inlet vel. 
(ft./sec) 


Coolant mass flow 
rata (107 lbse/hr) 


Head loss (ft H20) 
Pumping power (hp) 


Coolant temperature 
rise (9F) 


c ean temperature 
feeronce {OF} 


Heat transfer area 
(sq. ft.) 


Average Peg ett heat 
transfer co Sit 
ee aaaris 


xit- “fraction 
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Po TABLE V- 
Power Sinisization With Volume Constant 
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OPTIMIZED Signe { 
PARAMETER RESOLTS (%) 
Total number of tubes ! 5062 3.2 
| ® of tubes in cocler | 6.7 | -4.3 | 
Tube length (ft) 11.39 | +10.6 | 
Tube inner diam. (in) 594 } +12.7 
Tube wail thick. (in) 022 | -55.1 | 
Tube outer diam. (in) | -638 | 42.1 | 
Tube pitch 1.753 25.2 | 
| Void diameter (ft) ‘ 1.14 | 43.6 
Bundle diameter (ft) 6.73 | -6.1 | 
Condenser volume | 431.6 | -0.2 | 
(cu. ft.) 
i Area ratio | 0.297 #1127 | 
{ Coolant inlet vel. 6.74 -20.5 
(ft. /sec) 
Coolant mass flow 1.508 -2.1 | 
rate (107 lbm/hr) 
Head loss (ft H20) 4.69 ~36.2 | 
i Pumping power (hp) 34.77 ~37.6 | 
; | Coolant tenperature { 11.6 417.2 
: aie | ke | 
Lod mean temperature 27.63 -3.5 
ference (°F) | | 
1 Heat transfer area 9631.4 +9.3 
(sq.ft) | ( 
Average overall heat { 627.3 9.1 { 
transfer coefficient 
btu/ (hr-sq. ft.-9 PF) {| 
: es: exit- fraction 0.0 -100. | 
: (% of input) 
: Heat rejected 1.669 +15.0 
i (108 5 avees | 
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PARAMETER 


Total number of tubes 


% of tubes in cocler 
Tube length (f+) 

Tube inner diam. (in) 
Tube wall thick. (in) 
Tube outer diam. (in) 
Tube pitch 

void diameter (ft) 
Bundle diameter (ft) 


Condenser volume 
(cu. ft.) 


( 

| 

{ 

| 

Area ratios 
Coolant inlet vel. 
(ft. /sec) 

| Coolant mass fl 
rate (107 Tbayhr) 
Head loss (ft H20) 
Puaping power (hp) 


Coolan* temperature 
rise (°F) 


Loq sean temperature 
frerence (OF; 


Heat transfer area 
(sq. ft.) 


pver sds overall heat 
transfer coeff Gaene 
bta/ (hr-sq. f£t.- °F) 


+ ose exit-fracticn 


Ste 

(% of inpat) 
t 

A985 Beavaes¢ 


TABLE VI 
Volurae Minimization With Power Constant 


a cae PD SOP nt aA SP Py SD ct ON OS IED A eA a OD ae OI ST AE TS Ne OE OE I A SO ITT aaa 
- 
e 
wt 
= 
oe 
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Total number of tubes 
{ 


OE Se 


Voluse Minisgization With Linear Variations 


PARAMETER 


TABLE VII | 
{ 


% of tubes in cocler 


Tube length (ft) 


{ 
| | 
Tube inner diam. (in) | 2527 |+.725 583 | ---- | 
Tube wall thick. (in) 049 { 022 {755.1 | 
Tube outer diam. (in) 2625 |+.773 -627 | o--- 
Tabe pitch 1.4 *1.75 1.44 | ---- 
Void diameter (ft) 1.10 1.28 416.4 
Bundle diameter (ft) 7.17 6.8 -5.2 | 
Condenser volune 432.3 345.8 | -20.0 | 
(cu. ft.) ] 
Area ratios | 0.266 | 0.298 #12.0 | 
Coolant inlet vel. 8.473 9.0 46.2 
(ft. /sec} 
Coolant mass flow 1.540 | 2.48 +61.0 
rate (107 lbm/hr) ! 
Head loss (ft #20) | 7.35 5.84 | -20.5 { 
Pumping power (hp) } 55.69 71.1 427.7 
eoocens temperature 9.90 7.1 ~-28.8 | 
rise (°F) 
og mean tea perature 28.62 30.14 5.3 
ference (9%F) 
Heat transfer area 8814. 7746. -12.1 
(sg. ft.) 
Average overall heat 574.9 713.1 { 
transfer coefficient 
btu / (hr-sq. £t.-9F) 
Steam exit~-fraction 13.3 0.5 


(% of input) 
Heat rejected 1.451 1.665 
(108 bti/hr) 


* Inner row values followed by outer row values. 
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VI. CONCLUSIONS 


The intent of this research was to create a detailed 
condenser analysis ccde capable of being coupled with a 
numerical optimizer and to test the program to prove its 
versatility. An additional objective was +o validate the 
analysis with existing data. The results of the test cases 
ware presented in Chapter Pive; the rasulting conlusions are 
summarized here. 


A. There were significant difficulties encountered in 
formulating the complex condenser dasign anaysis, CONDIP, ina 
away that was compatible with the optinizer COPES/CONMIN. 
Howsver, the majority of those problems were cverconme 
resulting in the creation of a program which, when combined 
with the optimizer, is capable of taking any initial design, 
no matter how impractical or infeasible, and solving for an 
optinun solution based on a set of pre~determired 
constraints and design variables. There are still some 
minor limitations as to the degree of optimization, but the 
final design is usually within 10 pecent of the single best 
optimus. In additicn, the test cases indicate that as many 
as ten design variables and six constraints can be used 
Simultaneously in the design optimization with CONDIP. 


B. The test cases demonstrated the effectiveness of CONDIP 
as a design tool for not only the conceptual design of a 
condenser, but also in evaluating comparison studies based 
on any number of design variable coabinations. The nusber of 
possikle combinations of design objactives, design variables 
and design constraints implies liaitless possibilities to be 
explered and evaluated. 
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C. An attempt was made to verify CONDIP with existing data 
with inconclusive results. Part of the blame can be piaced 
on the rather inadequate quality and quantity of the data, 
but the general performance of CONDIP's condensar indicates 
a weakness in the analysis. As statad aarlier, the source of 
the this weakness may be found in the correlations used for 
condensate inundation. The constants used in the expression 
for correcting shell-side heat transfer coefficients are 
based somewhat on ccnjecture. Yet thay play a significant 
rola in the overall condenser performance. Despite this 
limitation, the ability to optimiza CONDIP*s detailed anal- 
ysis is a significant step forward over using the tradi- 
tional and limited HEI method. 


D. CONDIP incorporates features that further increase its 
appeal asa design tool. By possessing the ability to 
linearly vary pitch and tube diameter, a better under- 


standing of how to improve condenser performance based on 
its configuration is realized. The capability +o incorporate 
shell-side tube enhancement is another added plus. The 
possibilities that can now be investigated are limitless. 


| VII. RECOMSENDALIONS 


In addition to the insight that this investigation has 
given into the generation of automated condenser design 
5 programs, i+ has specifically addrassed the shortcomings and 
pitfalls which may be encountered along the way and offerad 

possible solutions to overcome tham. Presented below are 

recommendations for furthering the developmentof CONDIP as 

a completely versatile and accapted design pregram. 

A. Since tha weak link and the most signifcant unknown in 

condenser analysis is the effect of condensate rain in 
= typical condenser environments, subsequent research should 
be devoted to investigating this phenomenon and developing 
more precise analytic correlations. In particular, the 
effects of vwalocity and flow direction on she condensate 
film should be attended to. 

B. Perhaps in conjunction with the above, a test condenser 
should be constructed which is geometrically similar to the 
model proposed in CONDIP in order to physically observe and 
‘ record the condenser performance. This data could then be 
: used to either verify CONDIP or strengthen some of its 


anaysis. In addition, this condenser should be built such 
| that tha tube bank can be arranged in any combination of 
pitch, tube diameter and row spacing to fully appreciate the 
: effects of these variables. 
| C. A series of sensitivity studias should be conducted on 
CONDIP to fully exercise its capabilities and determine the 
telative effects of various designvariables on condenser 
performance. Tradeoff studies sisitlar to those performed in 
this research would be most beneficial to fully understand 
condenser behavior. 


Rae ee os es ne 


D. Additicnal subroutines should be created which would 
allow tube enhancement to be a design variable. This 


involves developing corralations between heat ‘zransfer 
enhancement and associated frictional losses. This type of 
relationship can be developed for both tube-side ard shell- 
side enhancement. 

BE. Finally, i+ is recommended that additional refinement 
be performed on the code to incraase its capability and 
flexibility. One such way is to somehow allow pitch, “ube 
diameter and tube wall thickness to vary linearly by row 
while still allewing the number of tubes to be a design 
variable. The options available are limitless. 


APPENDIX A 
GLOSSARY 


While it would mcst beneficial to present a complete 
glossary of all the variables used in CONDIP, the sheer 
number makes it difficult to presant a comprehensive list. 
However, CONDIP makes liberal use of comment cards to define 
aS many variables as possible to make the code e¢asier to 
follow. Therefore, the computer listing in Appendix cis 
available for reference. A list of the possible desiaqn vari- 
ables and constraints is provided here along with its corre- 
sponding position in the GLOBCS common block for easy 
reference in writing the appropriate COPES data cards. In 
addition, it will be spacified whether these variables can 
be used as design constraints or dasign variables. 


1. ALST: The length of the condenser and cooler tubes in 
feet. ALST is to be used only as a dasign variable. 


a 2. DELWP: The pressure difference between the inlet and 
* outlet coclant headers of the condanser bundle in psi. 
‘a DELWP is to be used only as a design variable. 


3. DELWPC: The pressure difference between the inlet and 
outlet coolant headers of the condenser bundle in psi. 
DELWPC is to be used only as a design constraint. 


4. GPLOW: The mass flow rate of the coolart in lbm/hour. 
GPLOW cannot be used as a design variable simultanouesly 
with VELBI. Otherwise it can be used as a design variable 
or a design constraint. 


5. SIDI: The tube inner diameter 9f the innermost rew of 


the condenser bundle in inches. SIDI is to be used cnly 45 
a design variable. 


| 7. SIDO: The tube inners diameter 9£f the outermost rew of 
the condenser bundle in inches. If there is to linear 
variation of tube inner diameter then this variable repre- 
sents the tube inner diameter of the entire ccndenser 
bundle. SIDO is to be used only as a design variable. 


8. PHP: The coolant pumping power in horsepower. PHP is to 


be2 used only as a design constraint. 


9. RSPA: The spacing between ccncentric crows in the con- 
denser burdle in inches. RSPA is to be used only as 2 
design variable. 


10. RADINS: The inner void radius of the condenser bundle 


in feet. RADINS is to be used only as a design variable. 


11. REWI: The tube-side Reynolds number of the coclant in 


the innermost row of the condenser bundle. REWI is to be 
used only as a design constraint. 


12. REWO: The tube-side Reynolds number of the coolant in 
the outermost row of the condenser bundle. If there is 29 


linear variation of tube inner diameter then this variable 
4 represents the tube-side Reynolds number of the entire 
condenser bundle. REWO is +o be used only as a constraint. 


* 13. SDDI: Tube pitch (tube spacing/tube outer diameter) of 
the innermost row cf the condenser bundle. SDDI is to be 
used only as a a design variable. 


14. SDDO: The tube pitch 9f the outermost row of the con- 

denser bundle. If there is no linear variation of tube 
pitch then this variable represents the tube pitch for the 
entire condenser bundle. SDDO is to be used only as a 
design variable. 


a “tera h Lesa Bhan 


15. SLDI: Ratio of tube length to tube outer diameter of 


the outermost row of «he condenser bundle. SLDI is to be 
used only as a design constraint. 


16. SLDO: Ratio of tube length to tube cuter diameter of 
the outermost row of the condenser bundle. If «here is no 
linear variation of tube pitch then this variable repre- 
sents the tube pitch for the entire condenser bundle. SLDO 


is to be used only as a design constraint. 


17. VELBI: The velocity of the coolant in feet/sec. VELBI 
cannet be used as a design variable simultanecusly with 
GFLOW. Otherwise it can be used as either a design 
constraint or as a design variable. 


18. XW1: The ratio of tube thickn2ss to tube inner diane- 
ter. XW1 can be used oily as a design variable. 


19. XW2: Tube thickness in inches. XW2 is to be used only 
as a design variable. XW2 and XW1 cannot be used sinultane- 
ously. 


20. VOL1: The overall condenser and cooler volume in cubic 


feet. VOL1 is to be used only as a d2sign constraint. 


21. VOL2: The volume occupied by the tube bank, excluding 
the volume of the inner void, in cubic feet. VOL2 is to be 
used only as a design constraint. 


22. INOTOT: The total number of tubes in the condenser and 
cooler combined. If Option 1 is being used then TNOTOT is 
+9 be used only as a design constraint. If Option 2 is 
being used then TNOTOT is to be used only as a design vari- 


able. 


23. BNDRAD: The condenser bundle in f2et. BNDRAD is ¢o be 
used only as a design constraint. 


24. ARATIO: The ratio of the total cross-sec*ional area of 
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the tubes (based on the tube outer diameter) +o the tube 
sheet area. ARATIO is to be used only as a design 
constraint. 


25. QODII: The tube cuter diameter of the innermost row in 
inches. ODII is to be used only as a design constraint. 


26. QDOZ: The tube cuter diameter of the outermost row in 
inches. If there is no linear variation of ‘tube inner 
G@iameter then this variable represants the ¢ube outer dias- 
eter of the entire condenser bundle. ODOT is to be used 
only as a design constraint. 


27. NVLECMAX: The maximum allowabla steam velocity into the 
cooler. VLCMHAX can be used only as a design constraint and 
only when a cooler is being designed in the systen. 


28. PRCCLR: The percent of the total number of tubes in the 
cooler. PRCCLR can be used only as a design variable. 


APPENDIX 8 
USERS MANUAL FOR CONDIP 


This appendix describes the data cards that are neces- 

a sary in order te couple any design prograa with 

COPES/CONMIN. Also described are cards illustrating data 
input required by CONDIP to initiate analysis. Thus, the 
data is divided into the COPES/CONMIN program section and 

the CONDIP-based condenser design program section. 
The COPES data is segmented into "blocks® for conven- 
* dence. All formats are alphanumeric for title, end and stop 
. . cards; F10 for real data; and I10 for integer data. The 
formatted input may be overridden by inserting ccasas 
between data entries. Comment cards aay be inserted anywhere 
. ia the data stack prior to the end card and are identified 
ii ty a dollar sign ($) in coluan 1. fhe COPES data stack aust 
terminate with an end card containing the word "END" in 
coluan 1-3. It should be noted that information pertaining 
only to single analysis and optimization is presented here. 
Information concerning the other options available in COPES 
along with further explanation of COIzS capabilities can be 

found in [Ref. 13]. 

The analysis data is also segmented into blocks for 
convenience and they begin immediately following the "“END®" 
card in the COPES @ata. No comment cards are pernitted here, 
and the analysis data stack aust tersinate with the word 
"STOP" in columns 1-4. This is where the initial design 
values are placed for entry into CONDIP. 

Default values are recommended for use in the following 
COPES data cards unless ctherwise noted. It is recossended 
that these values in the COPES data blocks be used until the 
user becoses fasiliar with the prograa. In addition a 


ahaa tae <1 eee : 
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sample data input is illustrated in figure B.1 at the end of 
this appendix. 


in? Ce VieRY nidh en I. SANRIO PSA 


PR SPT eR - 


| 
{ 


DATA BLOCK A 
DESCRIPTION: 
FORMAT: 20,6 


COPES Title Card 


: 1 2 3 4 5 6 7 8 
i eres {| } | 
| Cae a a ee ae ee ae ee 


BESABKS: 
1) This line is awailable for a briaf description. 


Le pres 


MOREA TION NE EIN Hs Cp RE IG EE 


DATA BLOCK 8 
DESCRIPTION: COPES Programs Control Parameters 
FOBMAZ: 7110 


1 2 3 4 5 6 7 8 


CONTENTS 


Calculation control 
yeaa. input and St?Ps Data of, block 
8 in 


red. Rea dat ° tIoA ° 
Onec cie through pel pcogeans Dhe Tonal 
Aas iis3 is reqaireds Beeatning data 


Ope liegt teh, pata of, Bloeks AE SE, 


tanger 2ariok ae opedaca sensitivity 


: Pield 1 determines program execution 


Tie as 3 are to be left blank for the CONDIP 
a cation “8. COPES /CONMIN. 


a ee ee ee ere eset eee Lee soe tenet ein tare oe meee 


DATA BLOCK ¢ 
DESCRIPTION: COPES Integer Optimization Control Parameters 


FORSAT: 8110 


1 2 3 4 5 6 8 
if lrean |urvoas HACHX1 | NFDG 


FIELD CONTENTS 

| 1 IPRINT: Broge control usei in optimization 
a 

| Q print uring optimization. 

{ 


ase in tia 1 and final optimization 
¢ abow u ction value and 
es: n v9. pEusf tics at each iteration. 
bove lus ‘eonstraine values, 
rect tee ee tor and move parameter at 


Hin if: ove fe fie if gradient Lafgc gat on. 
tor . ob je gob Propose constraints 


dria y erone; aynengione}. Sea op EfSnal. 


AEF TO EERE EeePC ACCC ST 
= Ww Ww 


| 2 ITMAX: Beyiaug 2 ber. of abst aization iterations 
i 3 ICMDI Rs: sepia: digegtion restart paraneter. 
| i poker sceting. “te 1oE cr of saagai cade ang. every 
NSCAL terati ch | - Segie es 
yarze dables accor ing to to oi f32: ng aalc s 
nput. DEPAULT = 
5 ITRA: Baaber of gabse ative iterations | which 
CORVergence cr rt H on Bistre. optfaization 
4 proces terminated. DEPAOLT 
iz © Maines EPR SEUSS SB GERI Ug isin +o 
* hth aser-*? ° the esign 
: DEbaur a. set fINGBs 
7 NACHX1: one the maxiagua n active 


| CB katiata ta ant ett oats a Gps = §DV+2. 
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DATA BLOCK ¢ (Continued) 


FIELD CONTENTS 
8 WPDG: Pinite difference gradient identifier. 
0 Hues inforaation is earn by 
1 * ob ject ive is comput 
3 kage BES of oh peers 
are coupy ted ay nit ai ft erence. 
2 gra $F formation is computed 
Tr tere 
BESARKS: 
y ff f WSCA oi ted @ ITRM = NACHX1 2 0 
| yatue of L= s suggested an 


HaTh TY. PATIT NT. PS, Seduced when the user becoms 


: The default vajues will be used if the card is either 
eft blan ue of zero is 


ora wa entered. 
oped eras’ t alge far infix. 2 ‘ EOD SP ae veiltr soe bd 
to nated prarature Recommended valué is ITMAX 


coaple it he condenser is ensure that no 
dhae 1698 ibe axe ie 
on, can nsi taseziy, fs ee ent on any coa bi- 
patio o3 foi -¢) var abies. Th s just s Pag the DEFAULT value 
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ay DATA BLOCK 2 


DESCRIPTION: COPES Floating Point Iptimization Progran 
Parameters 


FORMAT: 8F10 


FLELD CONTENTS 

| 1 FPDCH: Relative cpaave in design variables in 

\ Calcalating £3 finite difference gradients. 
EPAULT = 

2 PDCHN: Min yaa absolute, s in finite ds £ fference 

gra Tent ca i5n DEFAULT = 

3 CT: Constraint thickness parameter. 


eM ee POTTER AEE PO OT tae Ce 


| DEFAULT = -0.1. 
| 4 CTI: iY 


EFAULT = 5.0. 


inipua absolute value of CT considered 
nt ptiniz gation process. 
DEFAULT © 
ee | 5 CTL: ze onstygint +thicknass parameter for linear 
| e constraints. 
EE 6 CT LAIN: Piniaoa abso ute value of CTL considered 
opts Bization process. 
| } DEFAULT = 0.001 
7 THETA: Nea Bir e of h-off factor in th 
i § of feasibie Girest iors: . 
| § DEFAULT 
8 PHI: Participation coefficient ed if 
Bore ¢ chs fats are violated: lal 
i § 


t 
| 


DATA BLOCK D (continued) 
FORNAT: 2F10 


1 2 3 
DELFUN| DABFUN ee 


DELFUN: Minisum relative change in objective 
function o indicate convergence of 
optimization procass. DEFAULT = 0.001. 


DABFUN: inisga absojute change in objective. 
unction to icate Convergefice of the 
op tiaizat ion rocess. 
D AULF = 0-9 1 tiaes the initial 
objective value. 


RESARKES: 
1) Note that. dat xr Dat Block D i tered two 
separate cards. i f ank card 1eatatees the default value is 


f the NDV is ter t 3, ¢ lue f 
Bis eeeteeitaT5°tet ae Ns 22 Che seooesented falsentos 


DATA BLOCK & 


PTION: Total Number of Design Variables, Design 
BESCBL Objective Identification and Sign on Sasign 


Ob ject ive. 


2110, F10 


FORE AT: 


CONTENTS 


1 NDVTOT: Total number of variables linked to the 
: design variables. NDVTOT must be greater 


nore 
design varjable. the veige of eac 


pargaetar is the v e of the design 
variable times a aultiplier which may be 
{ G@ifferent for each parameter. 

| * DEFAULT = NDV. 

OBJ: Global va b aber associat with 

| ; : bg tive function 2 oper Trettoa — 

: Bum sensitivity analysis. 

GNOPT: Si ecti f izat 

2 oy toMhagnetey « eth tyeackion ts fa bee 

) max alze tT miniaized, +1.0 indicates 
maximization; -1.0 indicates ainimization. 


DEFAULT = -1.0 


BEBARKS: 


Currently there are not any variables in CONDIP which are 
Led fovaky gf the design ay arte? Therefore the DEFAULT 


ee Se et ae 


se a es Oe coat 


sear a 


ee ary eee 


RATA BLOCK £ 


DESCRIPTION: pesty3 Variable Bounds, Initial Values, and 
ag Pactors. 


Seal 


FORMAT: 4F10 


ps 
aa 


| 


FLELD CONTENTS 
VLE: Lower bound on the design variable. 
VOB: Upper bound on the design variable. 
3 : Initial value of the. des variable. 
r? x ts non-zero, this sift supercede the 
value initialized by subroutine ANALIZ. 
4 SCAL: Des wariable scale factor. Not 
bes§o? NSCAL 2 6 tn Bidck Cc 
RENABKS: 


here spust be one sepa e data card for each desi var- 
lhe @. therefore Bese Party be NDV data car — ony 


apolGE ar sdhlasGhiT'be "ei tered "e 
ne called in ANALIZ. 


initial values for the 
be entered through the INPUT subrou- 
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DATA BLOCK & 
DESCRIPTION: Design Variable Identification 


FORMAT: 2210,F10 


1 2 3 & 5 6 7 8 
osew | IDSGN | AMOLT 
| 
ey SEE | 
FIELD CONTENTS 

1 NDSGN: Design variable numbar associated 
with the variable. 

2 IDSGNs Global variable number associated 
with the variable. 

3 AMOLT: Constant multipliar on the variable. 
The value of 5 e variable will be the 
value of the sostan variable, NDSGN, tines 
AMOLT. DEFAULT = 

REMARKS: 


1) There aust be ong hae ayate card for each of the NDVTOT 

ables. es ata cards aust follow the same 
grces as - ao eveapondtes design variable parameter cards 
n oc ° 
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i 
+ RATA BLOCK # 

ag 4 DESCRIPTION: Number cf Constrained Parameters. 

- FORMAT: 110 

fo} 


= NCONS | 
4 __ se a A 


# FIELD CONTENES 
1 NCONS: Number of constraint SETS in the | 
optimization problean. P 
z REMARKS: 


If two or aore ad ent aramet2rs in the Global ccmmon 
block have the gage tL 2 P 
same@ constraint set. 


mposed, these are part of the 


DATA BLOCK 1 
DESCRIPTION: 
FORMAT: 3110 


1 2 3 


,- 1 ICCN: Pirst Global number corresponding to 
a the constraint set. 


2 JCON: Last Global nuaber correspondin 
the constraint sst. DEFAULT = I on? 


= 3 LCON: Linear congtraint i ngifier ger 
fy this set nett as variables. 
i LCCN = 1 ef dicates fed constraints. 
a DEFAULT = 0 2 Ronlizeas constraint. 


Constraint Identification and Bounds. 


CONTENTS 


RESARKS: 
2 Se COnDEE there if ek Global nuaber and thus one sd 

raint that couprise Seo straint set. Therefore @ 
DEPAULT value is used for J one 


2) All the constraints 4n this apalysis are nop} inear. Ths 
PAULT walue was therefore used- for LCON as wel 


3 ghis is the first card of a two card set which sust be 
read together. 


mpi Absa SA es dle 
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DATA BLOCK JX (Continued) 


FORMAT: 4F10 

1 2 3 4 5 6 7 8 
fea = 1 een | 
FLELD CONTENTS 

1 BL: Lower bound on tha constrained variables. 
Value less than -2.Q02+15 is assuned 
un bounded 

2 SCAL1: Normalization factor on lower b bound. 
DEFAULT = Max of ABS(BL) 

3 BO: Upper bound on tha const trajn ad variables. 
Value eater than +2.0E Ss assumed 
unbounded. 

4 SCAL2 No ization factor on uppe ound. 
DEPAUL? LT Max of ABS (BQ) BR 5.3° 

BESARKS: 
1) The able” azcepeicn tor ¢33 uspally be defa alt with 

e not cep ee 3f exit peccecn. normals Zation 

factor used his constraint usu aily® ten nes the 


upper bound. 


DATA BLOCK P& 
: DESCRIPTION: COPES Data ‘ZND* Card. 
| FORNAT: 341 


FIELD COUT ENTS 
1 The word *END* in coluan 1-3. 


1) This card gygt appear at the end of the COPES data. 
2) This ends the COPES input deck. 


eae we 
Lehi td 
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DATA BLOCK Aa 
DESCRIPTION: Geometry Option 
FORMAT: I5 


1 2 


nae egret te HEE ETT er OS 


CONTENTS 


condenser geometric options 

T = 1. Number of condenser rows is 

put as a constant; the number of 

Ss a_de Sndent variable. 

data bloc s EE,FF,and GG. 

T = Nuaber =f. tubes is allowed 

be an , independant variable and the 
rows is a dependant variable. 

eee b eck HH and I eb 

T va s IOPT = 2. 


OQp tHarte He 
wmNcIOONnNE Of 


BESARKS: 
2 Date is right-justified and blanks will be interpreted as 


ait efit data *3 ie ares tapit*, is because. tie ote = opt ion 


i; the des ign apairsis s less sensi the 
discussed earl Recon ended using the DEPAUL? va Lue of 
0.01 for FDCH. 


Rg RYE Bo 


Data BLOCK BB 
DESCRIPTION: Condenser Orientation 


FORNAZ: 125,3F10 


4 5 6 7 8 
veed_—f to 
ELELD CONTENTS 
1 ISEC The number of sectors in the condenser. 
2 SECWID Sector width in degrees of arc. 
3 Par Symusetry angle measure from tha vertical. 
ry PRCCLR The percent of tha tubes in the cooler. 


BESABRKS: 

iy De Data BB is required, no matter what geometry option is 
° lirsitation on ISEC and SECHID i at their - 

het ae 3 is than 36 Ce grees. product. S exact] P580 

tees iain nh trigonometric ace zons will retu 


Par is inet angle from the vertical that cuts the conden- 
zhr 


aL pete is right-justified and blanks will be interpreted as 
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DATA BLOCK Cc 
DESCRIPTION: Void Size, Tube Length, Row Spacing 
FORMAT: 3F10 


1 2 


RADINS 


cogent ADAP ERR AT EATON RES SRE RL” SIP 


CONTENTS 
RADINS The inner void radius; feet. 
ALST The tube length; feet. 
RSPA Goncentric tow spacing about the void; 
nche 


pie eipee? eens 


BESARKS: 


é <pate cc is required, no matter what geometry option is 


BATA BLOCK ODD 
DESCRIPTION: Tube Material Parameters 


s 15,4P10 


a Fig an dicecing the tube thickness 
ie icq ti 
Tube thickness is in 


tio of tube thiskness to tube’ lanes 


t 
aires 33. Tube thickness is input in 
nche 


The input $or wall thickness, dependent on 
the walue for IWALL. 


TOBESW $pesitc yeignt of the tube material; 
a/ (cu. 


SKW Tape 22te rial cheraa Bern eactee str: 


POOL Tube fouling factor. 


BESABKS: 


é pate DD is required, no matter what geometry option is 


2) Date is right-justified and blanks will be interpreted as 


Tae atta 


QATA BLOCK EE 
DESCRIPTION: Sumber of Rows 
FORSAT: 15 


FIELD CONTENTS 


1 NORQWS The nuaber of concentric rows in the 
congener bundle built around the center 


ERBABES: 
1) Data EE is used ggly vhen [opt = 4} 
2) Data is right-justified and blanks will be interpreted as 
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BEY Ye” ys 


min a pee 


DATA BLOCK FF 
DESCRIPTION: Tube Inner Diameter 


FORMAT: 115,2P10 


CONTENTS 


z1a9 to indicate whether tube irner 
ameter will linearly vary by row through 
e condenser un * 
diageter is uniform 
through the condenser bun 
pIAx z 2. Tube inner diageter ¢ varies 
aa ies through tha bundle by row. 


qube inner diameter of the outer row; 


ube inner diamatar of the inner TOW; 
nehes. 


BESARKS: 
1) Data PF is used only when IOPT = 1 


2 ba Data is right-justified and blanks will be interpreted as 


aceoces tubes use the inner diameter of the innermost bun- 


&) The DEFAULT value is NDIA = 1 


pare --3e_ ot 


DATA BLOCK && 
DESCRIPTION: Tube Pitch 
FORMAT: 115,2F10 
1 2 3 4 5 6 7 8 


MPITCH 


FIELD CONTENTS 
1 MPITCH =? ate to indicate whether tube pitch 
41} nearly. yary by row through the 
con en se 
1 MPITCH = 1, tape: pitch is uniform through 
the condenser bani} 
2 MPITCH = 2. Tube pitch _ waries linearly 
through the bundl3 by row. 
SDDO Tube pitch of the outer row; 
SDDI Tube pitch of tha inner row; 
BESARKS: 


1) Data GG is used only when IOPT = 1 
2 pata is right-justified and blanks will be interpreted as 


3) Ccoler tubes use the pitch of tha innersost bundle row. 
4) The DEPAUL? value is MSPITCH = 1 
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| 
: 


BATA BLOCK 8H 
DESCBEIPTION: Tube Inner Diametar and Tube Pitch 


FORMAT: 2F10 
3 4 5 6 7 


es ae a 
nas as ae fa 


CONTENTS 


1 SIDO Tube inner diameter for the entire 
condenser; inches 


on Tube pitch for tha entire condanser 


BESARKS: 
| 1) Data BH is used ogly when IOPT = 2. 
f 2) Data is right-justified and blanks will be interpreted as 


r 
In the calculations IDI is set equ t DO 
i ait i avoids the tased Jos tua eyacess 


set 9 ° oe f s 
nomenclature for each geometry option. 


BATA BLOCK IT 
DESCRIPTION: 
FORSAT: F2 


Total Nuaber of Tubes in the Condenser 


CONTENTS 


TNOTOT The total number sf tubes in the 


condenser (cooler and the sSundle). 


RESARKS: 


1) Data II is used only when IOPT = 2. 


a pete is right-justified and blanks will be interpreted as 


DATA BLOCK sy¥J 
DESCRIPTION: Inlet Steam Mixture. 


FORMAT: 15,2F10 
1 


lag indicating the type of 
“condensable gas entering the systen. 
S$ = 1. This indicates that the gas is 


¢ This indicates that the gas is 
n dioxide. . 
S = 3. This iniicates that the gas is 
ixture of the two. 


af flow rate entering the condenser; 


° 


atjio_of the non-sond 
o inlet steam flow; 


ensable gas flow 
lbm/hr. 


BESARKS: 


1 guess JJ is required, no matter what geometry option is 
chosen. 


a) pak is right-justified and blanks will be interpreted as 
e * 


DATA BLOCK kK 
DESCRIPTION: Inlet Temperatures. 
FORMAT: 2F10 


1 2 3 4 
Pati ile 
szez_{ssum 


STBI Coclant inlet temperature; OF. 


CONTENTS 


STSAT1 Talet steam saturation temperature; 


REMARKS: 


3 seek KK is reguired no matter, what geometry option is 
en. 


2) 
ze 


Date is right-justified and blanks will be interpreted as 


rt 


a et 


a erre—arrers 


RATA BLOCK LL 


DESCRIPTION: Cooling Water Parametars 


FORMAT: 2F 10 
teow | 
FLELD CONTENTS 
1 IFLOW A ccntrol flag for cooling water 
specifications. 
1 IPLOW = 1. Input pressure: ¢ aFoP across 
cooling water heaiers in ia. 
2 IPLOW = 2. Input cooling pater velocity 
in Eis secs 
3 IFL 3. Input coolant flow in lbm/hr. 
2 x5 Actually input the value for flicw into 
fis vari able. specctication fo for flow 
to he ieterained ar 
REMARKS: 
bate LL is required, no matter what geometry option is 
sen. 


a pare is right~justified and blanks will be interpreted as 


3) X5 acts as a temporary all-purpose storage variable for 
atever expression for coolant flow is used. 


I AS 86 oe a Ok OE Ele NN = 
: ar “AY ye EW oe ee 


« 


DATA BLOCK &8 


DESCRIPTION: Internal Enhancement Ragions 


FORMAT: 15 
1 2 3 4 5 6 7 8 
=| | [_ ees 
FIELD CONTENTS 
1 NEI Nusber of internal enhancement regions. 
A walue between 1 ani 6. 
REMARKS: 


1) Data MM is optional. There must be NEI Subs equept data 
cards providing the necessary paramaters for each tegion. 


2) enhancement can only be used if [opT = 1. 
3) Data is cight-justified and blanks will be interpreted as 


&) This value was zero for all runs. 
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DATA BLOCE Sy 
DESCRIPTION: Internal Enhancement Parameters 


FORMAT: 225,3P10 


8 
ELELD CONTENTS 

1 NRNI Row nusber of first row in internal 
enhancerent region. 

2 NETI Nuaber of tubes in each internal 
enhancement region. 

3 ENI Internal heat transfer enhancement 
factor. 

4 BC Coefticie nt in iaternall enhpansed “ube 
coolant pressure irop calcu ation. 

5 BE Expenent in coolant pressure dro 
calcataticn. P P 


REMARKS: 

NF Data NN is optional. Hovever, it NEI i 
then there must NEI "NN" data cards 

data for each enhancesaent region. 
2) enhancement can only be used if IOPT 2 1. 


a Data is right-justified and blanks will be interpreted as 


4) Thege values are Gonetent for entire run and cannot he 
changed by the optinizer. 


5) This value was zere for all runs. 


s 9 ter than zero 
to provide the necessary 
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DATA BLOCK 09 
DESCBIPTION: External Enhancemant Regions 
FORMAT: I5 

1 


CONTENTS 


Number of external enhancemant regiors. 
A value between 1 and 6. 


BESARKS: 


1) Data oOo is optional. There aust be NgE subsequent data 
cards providing the necessary parameters for each region. 


2 Data is right-justified and blanks will be interpreted as 


4) This value was zere for all runs. 


a mrp MOR We Swe 
re re eat 


Rees 


Niece, is agtaesioxcnen : 2 . SAN mars Ae Na et So 


DATA BLOCK PP 
DESCRIPTION: External Enhancesen* Parameters 


FORMAT: 225,2F10 


5 6 q 8 
{| 

FLELD CONTENTS 

1 WRNE Row nuaber of first row in external 

P enbancesent region. 

2 NETE Number of tubes in aach external 

enhancement region. 
3 ENT yx erat heat transfer enhancement 
actor. 
Hy ENH Stean-side pressure drop factors. 


eae 
Data PP is optiona verer iF NEE is eater than zero 
then Rete aust HEE Med" data csrds provide the necessary 
ta for each én hancenent Fegion. 


2) Data is right-justified and blanks will be interpreted as 


3) Enhancement can only be used if I0PT = 1. 


ase? values are. enstant for an entire run and cannot be 


chan by the opt zer. 


5) This walue was zero for all runs. 


a cece, 


DATA BLOCK 00 
DESCRIPTION: Baffle Options 
FORSAT: 15 


FIELD CONTENTS 
1 IBAF A flag to be used to determine baffle 


nusaber and location. 


BESARKS: 
1) Data QQ is optionaL 
2) bata is right-justified and blanks will be interpreted as 


3) These values are constant for an entire run and cannot be 
changed by the optimizer. 


4) adds tjonal speci tied baffles were not used in any of the 
tuns. This value was zero for all runs. 
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DATA BLOCK &8 
DESCRIPTION: Baffle Location 


FORMAT: 15 
1 2 3 4 5 6 7 8 
JBAF | | 
FIELD CONTENTS 
1 JBAP An array containing baffle locations 
BESARKS: 


1) Data RR is optional. 

2) bata is right-justified and blanks will be interpreted as 
3} Thege valves ar? gonstant for an entire run and cannot be 
changed by ¢ optisizer. 


4) Additional specified baffles were not used in any of the 
runs. This value was zero for all runs. 
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DATA BLOCK $8 
DESCRIPTION: 


FORSAT: 15 
1 2 3 


Detailed Printout 


4 5 6 7 8 
Io 


CONTENTS 


1 IPRT A flag to generate a detailed output 
of the condenser analysis (OUT3) 


a 


RBEGARKS: 


1) Data SS is optional. 

" 2) Data is right-fustified and blanks will be interpreted as 
a 3} Thege valugs are gonscant for an entire run and cannot be 

i changed by the optisizer. 


4 This value was zero for all runs. 


e ‘. 
a Ne ee See te 
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APPEBDIZ C 
CONDIP LISTING 


The following Appendix contains a complete listing for 
CONDIP. An effort has been made to aake the program as reaad- 
ables as possible through liberal usa of comment cards. 
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